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Abstract 


Non-conventional  bio-derived  fuels  have  been  evaluated  for  use  in  hybrid  rocket 
motors.  Tests  were  conducted  at  combustion  pressures  in  the  range  of  100  -  220  psig  and 
thrust  levels  of  40  -  170  newtons.  Beeswax  was  tested  with  oxygen  as  the  oxidizer  and 
showed  a  regression  rate  at  least  three  times  as  high  as  traditional  hybrid  propellant 
combinations  such  as  hydroxyl-terminated  polybutadiene  (HTPB)  and  liquid  oxygen 
(LOX).  This  provides  the  promise  of  a  high  thrust  hybrid  rocket  motor  using  a  simple, 
single  port  geometry  and  overcomes  the  main  weakness  of  traditional  hybrid  rocket 
motor  propellants,  which  are  low  regression  rates.  Beeswax  was  also  tested  with  nitrous 
oxide  as  an  oxidizer,  but  further  testing  is  needed  to  attain  high  enough  combustion 
chamber  pressures  to  achieve  stable  combustion.  Experimental  evaluation  of  the  specific 
impulse  for  beeswax  and  oxygen  was  moderately  successful  for  lab  scale  testing,  but 
needs  further  refinement.  Analytical  studies  were  performed  to  evaluate  the  theoretical 
performance  of  non-conventional  hybrid  rocket  motors.  This  analysis  indicates  beeswax, 
lard,  a  mixture  of  paraffin  and  lard,  and  combinations  of  beeswax  and  aluminum  should 
all  perform  better  than  traditional  hybrid  rocket  propellants  considered  when  burned  with 
oxygen.  For  a  combustion  chamber  pressure  of  500.38  psig,  beeswax  and  oxygen  yielded 
a  maximum  specific  impulse  of  327  s.  The  high  specific  impulse  combined  with  a  high 
regression  rate  combine  to  make  beeswax  and  oxygen  a  potentially  high  performing 
hybrid  rocket  motor  propellant  for  launch  vehicles,  suborbital  rockets,  or  orbital  kick 
motors. 
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I.  INTRODUCTION  AND  BACKGROUND 


I.I  Background  on  Hybrid  Rockets 

Hybrid  rocket  motor  technology  has  been  around  since  the  1930s.  A  hybrid 
rocket  motor  typically  uses  a  solid  fuel  and  a  liquid,  or  gas,  oxidizer.  There  are  several 
reasons  hybrid  rocket  motors  are  desirable  compared  to  solid  rocket  motors  or  liquid 
rocket  engines.  Better  performance  can  sometimes  be  achieved  by  propellants  that  are  in 
different  states  (i.e.  liquid/solid  or  gas/solid).  Some  solid  fuels  that  can  be  used  in  hybrid 
rocket  motors  can’t  be  used  in  solid  rocket  motors  due  to  incompatibility  with  other  solid 
fuels  and  oxidizers.  Hybrid  rocket  motors  are  safer  than  solid  rocket  motors  because  it  is 
virtually  impossible  to  have  a  combination  of  the  fuel  and  oxidizer  that  can  result  in  an 
explosion  due  to  the  unique  combustion  process  in  a  hybrid  motor.  If  a  crack  develops  in 
the  solid  fuel  grain,  the  combustion  process  will  not  allow  a  dramatic  increase  in  burn 
surface  area  as  occurs  in  a  solid  rocket,  which  in  turn  prevents  a  continuously  increasing 
chamber  pressure  that  would  eventually  result  in  an  explosion.  This  is  due  to  the 
combustion  occurring  in  a  diffusion  flame,  as  shown  in  Figure  1. 
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Figure  1 :  Hybrid  Rocket  Motor  Combustion  Process 
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The  fuel  enters  the  boundary  layer  by  vaporization  at  the  solid  fuel  surface,  and  the 
oxidizer  enters  the  boundary  layer  from  the  free-stream  flow  by  diffusion.  When  the  fuel 
and  oxidizer  meet  at  approximately  the  stoichiometric  ratio,  a  combustion  zone  occurs. 
Due  to  the  diffusion  flame  occurring  in  the  boundary  layer,  there  is  no  increase  in  surface 
area  when  a  crack  forms  in  the  solid  fuel  grain  because  there  is  little  to  no  oxidizer  within 
the  crack.  Finally,  unlike  solid  rocket  motors,  a  hybrid  can  be  shutdown,  restarted,  and 
throttled.  In  general,  traditional  choices  for  propellants  in  hybrid  rockets  have  one 
primary  weakness.  Due  to  the  regression  rate  being  low  as  compared  to  solid  rocket 
propellants,  a  high  thrust  hybrid  motor  requires  multiple  combustion  ports  to  produce 
thrust  levels  similar  to  those  of  a  solid  rocket  motor  using  a  single  combustion  port.  This 
is  necessary  because  the  thrust  produced  by  the  motor  is  a  function  of  the  propellant  mass 
flow  rate,  which  is  a  function  of  the  fuel  regression  rate.  With  a  low  regression  rate,  a 
motor  needs  more  combustion  ports  to  increase  the  propellant  mass  flow  rate  to  produce  a 
high  thrust  level.  This  results  in  increased  complexity  in  the  solid  grain  fabrication  and 
poor  volumetric  efficiency.  Unlike  liquid  rocket  engines,  the  fuel  mass  flow  rate  in  a 
hybrid  rocket  motor  depends  on  the  oxidizer  mass  flow  rate  through  the  combustion  port. 
The  thrust  of  the  motor  also  depends  on  the  specific  impulse,  which  is  a  function  of  the 
combustion  chamber  temperature.  Another  negative  aspect  of  hybrid  rocket  motors  is 
their  unproven  flight  heritage  compared  to  liquid  and  solid  rockets  (1,2). 

I.II  Review  of  Outside  Hybrid  Rocket  Research 

Over  the  years,  there  has  been  an  extensive  amount  of  work  on  hybrid  rocket 
motor  technology.  A  significant  amount  of  this  work  has  been  on  improving  the 
performance  of  hybrid  rockets  in  order  to  develop  an  efficient,  high  thrust  rocket  motor 
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with  a  simple  grain  geometry.  Hybrid  rockets  have  been  considered  for  uses  in  launch 
vehicles,  sounding  rockets,  and  satellite  propulsion  systems.  In  fact,  due  to  the  extremely 
safe  nature  of  hybrid  rockets,  they  are  being  used  extensively  in  the  amateur  rocket 
community. 

One  typical  goal  of  hybrid  rocket  research  is  to  evaluate  the  regression  rate  as  a 
function  of  the  oxidizer  mass  flux  rate  (1).  The  oxidizer  mass  flux  rate  is  defined  as  the 
ratio  of  the  oxidizer  mass  flow  rate  to  the  combustion  port  cross-sectional  area  as  shown 
in  Equation  1 

m 

(1) 

Ap 

where  Gox  is  the  oxidizer  mass  flux  rate,  mox  is  the  oxidizer  mass  flow  rate,  and  Ap  is  the 

combustion  port  cross-sectional  area.  The  regression  rate  can  be  expressed  as  a  function 
of  the  oxidizer  mass  flux  rate  as  shown  in  Equation  2 

r  =  aGj'  (2) 

where  r  is  the  regression  rate  and  a  and  n  are  experimentally  derived  constants. 

Research  at  the  University  of  Surrey  has  shown  that  hybrid  rockets  can  be  a 

feasible  option  for  satellite  propulsion  (3).  A  satellite  requires  a  propulsion  system  for 

orbit  maneuvering,  orbit  maintenance,  and  stationkeeping  maneuvers.  It  has  been 

demonstrated  that  a  hybrid  rocket  system  using  85%  hydrogen  peroxide  (HTP)  as  the 

oxidizer  and  polyethylene  as  the  solid  fuel  can  serve  as  a  cost-effective  solution  to 

satellite  propulsion  requirements.  The  hydrogen  peroxide  was  passed  through  a  catalyst 

to  decompose  it  into  steam  and  oxygen  at  a  temperature  of  approximately  630  °C  before 

injection  into  the  combustion  port.  This  hybrid  rocket  motor  displayed  regression  rates 
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similar  to  hydroxyl-terminated  polybutadiene  (HTPB)  and  liquid  oxygen  (LOX). 
Combustion  chamber  pressures  achieved  were  250  psi.  It  was  shown  that  a  hybrid  rocket 
motor  can  be  a  safe,  reliable,  and  cost-effective  alternative  to  a  solid  rocket  motor  for  an 
upper  stage  to  boost  a  satellite  into  a  higher  orbit.  Work  at  Purdue  University  has  also 
shown  that  80-88%  hydrogen  peroxide  and  polyethylene  can  function  as  an  effective 
hybrid  rocket  system.  The  tests  at  Purdue  were  conducted  over  a  range  of  combustion 
chamber  pressures  of  100  -  400  psi,  but  a  typical  test  chamber  pressure  was  200  psi.  As 
with  the  tests  at  Surrey,  they  used  a  catalyst  pack  to  decompose  the  HTP  for  the  ignition. 
This  type  of  process  provides  a  self-ignition  behavior  that  precludes  the  need  for  an 
ignition  system  and  provides  a  reliable  and  repeatable  ignition  process.  Smooth,  stable 
combustion  with  90%  efficiency  has  been  achieved  (4).  Amateur  rocket  enthusiasts  have 
focused  on  hybrid  motors  for  use  in  suborbital  rockets  all  over  the  world  due  to  their 
inherent  safety.  The  MARS  Amateur  Rocket  Society  (MARS)  in  England  has  developed 
hybrid  rocket  motors  for  sounding  rockets  using  nitrous  oxide  as  the  oxidizer  and 
polyethylene  as  the  fuel  (5).  Their  B4  hybrid  motor  has  been  successfully  static  tested 
and  produced  a  thrust  level  of  3,000  newtons  (N).  They  have  tested  a  down  rated  version 
of  the  B4  (1,400  N),  and  it  achieved  an  altitude  of  approximately  one  mile.  The  B4 
hybrid  rocket  motor  is  being  designed  as  an  upper  stage  for  MARS’  Deimos-3  hybrid 
rocket.  The  Deimos-3  is  being  designed  as  a  two  stage,  hybrid  rocket  that  will  be  able  to 
put  a  small  payload  to  an  altitude  of  62  miles,  which  is  the  edge  of  space.  The 
combination  of  polyethylene  and  nitrous  oxide  exhibits  the  same  low  regression  rate  as 
HTP  and  polyethylene,  so  it  is  difficult  to  obtain  high  thrust  levels.  MARS  has  the 
distinction  of  launching  the  first  hybrid  rocket,  amateur  or  professional,  from  British  soil. 
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Orion  Propulsion  has  produced  a  hybrid  motor  for  the  X-Prize  competition  that  uses 
asphalt  as  the  solid  fuel  and  nitrous  oxide  as  the  oxidizer  (6).  The  hybrid  motor  was  be 
mounted  on  a  Chevy  pickup  truck  for  demonstration  purposes,  and  has  a  maximum  thrust 
of  3,000  pounds  force  (lbf)  with  a  maximum  combustion  chamber  pressure  of 
approximately  400  psig.  Using  this  hybrid  rocket  motor,  the  truck  has  been  accelerated 
from  0  to  60  miles  per  hour  in  4. 18  s.  The  hybrid  motor  has  an  oxidizer  to  fuel  ratio  of 
6.0,  an  average  thrust  level  of  1,725  lbf,  an  average  combustion  chamber  pressure  of  207 
psig,  and  a  burn  time  of  10  seconds.  One  method  attempted  to  increase  the  regression 
rate  in  hybrid  rocket  motors  has  been  to  introduce  aluminum  into  the  solid  fuel.  At  Penn 
State  University,  nano-sized  particles  of  aluminum  powder  with  diameters  of  70  -  150 
pm  were  mixed  into  HTPB  fuel  grains  and  tested  (7).  Pure  oxygen  was  used  as  the 
oxidizer  at  combustion  chamber  pressures  up  to  1,750  psig.  Typical  tests  were  run  at 
oxidizer  mass  flux  rates  of  140  -  850  kg/m2s  and  chamber  pressures  of  320  -  650  psig. 
The  result  of  these  tests  showed  the  regression  rate  increased  by  50%  over  the  HTPB 
without  aluminum.  At  Stanford  University,  research  has  been  conducted  on  solid  fuels  in 
the  hope  of  producing  higher  regression  rates  than  previous  hybrid  rocket  fuels,  such  as 
HTPB  and  Plexiglas  (8).  This  was  based  on  previous  research  by  the  U.S.  Air  Force  into 
solid  cryogenic  fuels  such  as  frozen  pentane  and  frozen  oxygen.  Pentane  is  an  alkane 
hydrocarbon  that  is  used  as  a  fuel  or  a  lubricant.  These  tests  showed  frozen  pentane 
exhibited  regression  rates  three  to  four  times  higher  than  with  HTPB.  The  hypothesis  for 
the  increased  regression  rates  was  a  lower  heat  of  vaporization  of  the  cryogenic  solid 
fuels.  Stanford  developed  a  liquid  layer  hybrid  combustion  theory  in  an  attempt  to 
explain  the  higher  regression  rates  using  frozen  pentane  (9,10).  It  was  determined  that 
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the  liquefying  nature  of  the  fuel  was  the  primary  characteristic  that  contributed  to  the 
increased  regression  rates.  This  led  to  the  choice  of  paraffin  as  the  solid  fuel,  and  when 
combined  with  gaseous  oxygen  as  the  oxidizer,  yielded  significantly  higher  regression 
rates  than  HTPB  in  lab  tests.  In  larger  scale  tests,  it  has  been  used  in  a  10  foot  long 
sounding  rocket  launched  to  an  altitude  of  16,000  feet  in  Nevada.  They  have  conducted 
over  250  lab  tests  to  evaluate  the  regression  rate  of  paraffin  and  test  fired  motors  with 
thrust  levels  up  to  2,500  N.  The  proposed  mechanism  behind  these  higher  regression 
rates  is  the  formation  of  a  melting  layer  at  the  burning  surface  (11-14).  This  melt  layer  is 
hydrodynamically  unstable,  which  results  in  droplets  of  the  fuel  being  entrained  in  the 
combustion  zone  as  well  as  the  normal  gasification  of  the  solid  fuel.  One  conclusion 
from  the  research  was  that  the  surface  tension  and  melt  layer  viscosity  play  an  important 
part  in  determining  how  well  the  entrainment  of  the  liquid  fuel  droplets  occurs.  For  high 
regression  rates,  the  desired  characteristics  are  low  surface  tension  and  low  viscosity. 
Fuels  like  HTPB,  polyethylene,  and  Plexiglas  have  high  surface  tension  and  viscosity, 
which  explains  why  they  have  lower  regression  rates  than  paraffin.  Paraffin  and  oxygen 
yielded  regression  rates  3  to  4  times  higher  than  HTPB.  Further  research  on  paraffin  as  a 
high  regression  rate  hybrid  fuel  has  been  conducted  at  the  United  States  Air  Force 
Academy  (15).  Paraffin  was  tested  with  90%  pure  hydrogen  peroxide  as  the  oxidizer. 
Small  scale  lab  tests  were  conducted  with  the  hydrogen  peroxide  flowing  through  a 
catalyst  to  decompose  it  into  steam  and  oxygen  at  a  temperature  of  800  °C  where  it  auto- 
ignited  the  fuel.  The  tests  resulted  in  regression  rates  of  2.87  -  5.28  mm/s  for  oxidizer 
mass  flux  rates  of  1 1 1  -  162  kg/nrs.  These  regression  rates  were  higher  than  with 
HTPB/LOX  and  even  paraffin  and  oxygen.  For  a  oxidizer  mass  flux  rate  of  130  kg/m  s, 
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paraffin  with  hydrogen  peroxide  had  a  regression  rate  of  3.23  min/s  where  paraffin  with 
oxygen  yielded  a  regression  rate  of  only  2.6  min/s.  It  is  theorized  that  the  increased 
temperature  of  the  oxygen  after  the  hydrogen  peroxide  decomposes  is  the  main  driver  of 
the  increase  in  the  regression  rate  over  using  paraffin  with  gaseous  oxygen.  This  would 
indicate  that  other  materials  with  similar  physical  properties  as  paraffin  could  also 
potentially  have  equally  as  high  regression  rates. 

I  .III  Review  of  University  of  Tennessee  Hybrid  Rocket  Research 

Work  on  hybrid  rockets  at  the  University  of  Tennessee,  Knoxville  (UTK)  began 
in  the  fall  of  1999  (16).  A  project  was  funded  by  Lockheed  Martin  to  develop  a  hybrid 
rocket  test  stand  and  initial  research  into  hybrid  rocket  propulsion.  A  test  stand  was 
constructed  by  a  group  of  undergraduate  students  working  under  Dr.  Evans  Lyne,  and 
tests  were  run  using  polymethyl  methacrylate  (PMMA)  as  the  solid  fuel  and  gaseous 
oxygen  as  the  oxidizer.  It  utilized  an  ignition  system  using  an  injection  of  propane  in 
conjunction  with  the  firing  of  a  spark  plug  to  provide  the  initial  heat  of  vaporization  to 
the  solid  fuel.  The  test  stand  was  used  successfully  to  test  fire  hybrid  rocket  motors  with 
combustion  chamber  pressures  ranging  between  50-150  psig. 

Further  work  at  UTK  progressed  into  working  with  paraffin  and  lard  as  the  solid 
fuels  (17).  Following  on  the  work  at  Stanford  University,  lab-scale  tests  were  conducted 
to  verify  the  regression  rates  reported  for  paraffin  and  oxygen.  These  tests  were 
successful  in  showing  the  elevated  regression  rates  previously  documented.  Testing  then 
proceeded  with  lard  to  determine  its  suitability  as  a  hybrid  fuel.  Compared  to  paraffin, 
lard  has  a  relatively  low  melting  temperature,  and  this  caused  some  problems  with  its 
physical  stability.  During  testing,  due  to  the  fairly  low  combustion  pressures,  spallation 
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of  unburned  lard  through  the  nozzle  occurred.  This  led  to  mixing  the  lard  with  paraffin  in 
an  attempt  to  improve  the  mechanical  stability  of  the  solid  fuel  grain.  It  was  hoped  that 
the  increased  melting  temperature  of  a  50%/50%  mixture  by  weight  would  improve  the 
mechanical  properties  of  the  fuel  grain  and  reduce  spallation.  Despite  this  mixture,  the 
fuel  had  to  be  stored  in  a  freezer  until  the  time  of  firing.  Results  of  the  testing  showed  the 
lard/paraffin  mixture,  when  burned  with  oxygen,  showed  a  significant  increase  in 
regression  rate  over  HTPB/oxygen.  Another  attempt  to  improve  the  mechanical  stability 
of  lard  used  a  sponge  soaked  in  lard  as  the  solid  fuel  grain.  Once  the  sponge  was  fully 
saturated  with  lard,  it  was  placed  in  a  freezer  to  solidify.  This  seemed  to  work  fairly  well 
in  keeping  the  lard  from  leaving  the  nozzle  without  having  been  burned. 

I.IV  Motivation  For  Current  Research 

There  are  several  motivations  to  pursue  a  non-toxic  propellant  combination  for 
use  in  hybrid  rocket  motors.  One  is  to  find  a  propellant  combination  that  has  less  toxic 
exhaust  products  when  compared  to  some  solid  and  liquid  rockets.  For  example,  a 
typical  Ariane  5  launch  releases  270  tons  of  concentrated  hydrochloric  acid  into  the 
atmosphere  per  launch  (18).  Rockets  using  unsymmetrical  dimethyl  hydrazine  (UDMH) 
release  so-called  nitrogen  radicals,  which  are  potentially  damaging  to  the  ozone  layer. 
Also,  some  oxidizers  (i.e.  HNO3,  F2,  C1F3)  are  toxic  to  handle,  and  therefore  require 
additional  safety  precautions  during  fabrication  and  loading,  which  increases  the  cost  of 
hybrid  motor  production.  Another  consideration  is  the  contamination  of  groundwater  due 
to  the  chemicals  used  in  rocket  propellants  (19,20).  Perchlorate  is  commonly  used  in 
solid  rocket  motors  in  the  form  of  ammonium  perchlorate.  Perchlorate  has  been  found  to 
be  linked  to  thyroid  problems,  especially  in  children.  Many  military  sites  have  a 
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significant  amount  of  perchlorate  contamination  and  the  cleanup  will  come  at  significant 
financial  expense.  Civilian  rocket  manufacturers  have  also  felt  the  cost  of  cleaning  up 
perchlorate  contamination  (21).  Aerojet  has  spent  an  estimated  $250  million  over  the 
past  25  years  to  reduce  the  levels  of  perchlorate  contamination  to  groundwater  at  its 
Sacramento  facility.  Aerojet  also  had  to  pay  a  $25  million  settlement  due  to  a  jury 
decision  finding  them  responsible  for  the  deaths  of  three  individuals  and  the  illnesses  of 
four  others  due  to  drinking  tap  water  contaminated  with  perchlorate.  The  three  deaths 
were  the  result  of  lymphoma  and  the  four  who  survived  had  a  thyroid  disease.  Studies 
have  also  shown  that  Lockheed  Martin  is  responsible  for  perchlorate  contamination  to 
lettuce  grown  using  water  from  the  Colorado  River,  which  has  been  contaminated  over 
the  years  by  Lockheed  Martin  (22).  This  could  affect  what  the  FDA  decides  as  to  what 
level  to  set  for  safe  exposure  to  perchlorate.  Knowing  that  tap  water  and  lettuce,  along 
with  other  vegetables,  exposes  people  to  perchlorate  could  lower  the  threshold  they  set. 
The  combination  of  environmental  and  safety  issues  emphasize  the  need  for  a  propellant 
that  is  non-toxic,  and  has  relatively  benign  combustion  products.  At  the  same  time,  the 
motor  should  achieve  a  high  regression  rate  and  a  high  thrust  level  utilizing  a  simple, 
single  combustion  port  design  similar  to  solid  rocket  motors. 

Beeswax  was  singled  out  as  a  candidate  for  our  solid  fuel  for  several  reasons. 
First,  it  is  very  similar  to  paraffin  and  could  potentially  have  a  comparable  regression 
rate.  Second,  as  a  natural  substance  produced  by  bees,  it  is  non-toxic  and  can  be  easily 
obtained  in  mass  quantities  from  commercial  companies.  Third,  it  is  inexpensive,  widely 
available,  and  easily  formed  into  fuel  grains.  A  hybrid  rocket  system  using  beeswax  as 
the  solid  fuel  would  have  a  lower  cost  compared  to  other  rocket  systems  since  there 
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would  be  no  requirement  for  elaborate  safety  precautions.  When  burned  with  oxygen  or 
nitrous  oxide,  some  of  the  combustion  products  are  carbon  dioxide,  water,  carbon 
monoxide,  nitrogen  oxide,  etc.  These  exhaust  products  are  similar  to  the  emissions  from 
automobiles,  which  are  considerably  less  toxic  than  other  substances  released  into  the 
atmosphere  from  some  rockets. 

Two  oxidizers  were  considered  for  this  research.  Gaseous  oxygen  was  chosen 
because  it  is  easily  obtainable,  has  a  history  of  good  performance  as  an  oxidizer,  and  is 
easy  to  handle.  Also,  based  on  the  research  at  Stanford  with  paraffin  and  gaseous 
oxygen,  it  seemed  a  reasonable  choice  for  use  with  beeswax.  Nitrous  oxide  was  also 
considered  because  it  is  easy  to  obtain,  safe,  has  so-called  “self-pressurizing” 
characteristics,  and  is  commonly  used  in  hybrid  rockets  around  the  world. 

This  present  research  will  have  both  an  analytical  and  an  experimental 
component.  There  are  four  major  parameters  to  be  determined  in  the  experimental 
research  program.  The  most  important  goal  is  to  determine  the  regression  rate  of  non- 
conventional,  bio-derived  propellants  and  determine  functional  relationships  between  this 
rate  and  various  test  conditions.  An  understanding  of  the  regression  rate  for  a  fuel  is 
critical  when  designing  a  hybrid  rocket  motor  to  produce  a  desired  thrust  level. 

Additional  goals  will  be  to  measure  combustion  chamber  pressure,  thrust  levels,  and  flow 
rates  for  the  fuel  and  oxidizer.  These  measurements  will  allow  the  calculation  of  specific 
impulse.  The  main  objectives  of  the  analytical  work  are  to  obtain  theoretical  values  of 
specific  impulse,  thrust,  and  other  performance  parameters,  compare  them  to  the 
experimental  values,  and  determine  the  theoretical  performance  of  a  wide  variety  of 
propellants  over  a  range  of  equivalence  ratios  to  help  explain  the  experimental  results. 
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The  result  of  this  analytical  work  will  be  the  formation  of  a  repository  of  thermodynamic 
and  thermophysical  properties  for  previously  not  considered  fuel/oxidizer  combinations 
that  may  be  used  by  future  designers  of  hybrid  rocket  motors. 

The  analytical  approach  and  methods  will  be  outlined  in  Chapter  2.  The 
analytical  research  results  will  be  discussed  in  Chapter  3.  The  experimental  setup  and 
procedures  for  the  experimental  testing  will  be  discussed  in  Chapter  4.  Evaluation  of  the 
results  of  the  experimental  research  effort  and  a  comparison  to  the  some  of  the  analytical 
results  will  be  presented  in  Chapter  5.  Some  conclusions  and  recommendations  for  future 
work  will  be  presented  in  Chapter  6 
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II.  DESCRIPTION  OF  ANALYTICAL  RESEARCH  PROGRAM 


The  analytical  portion  of  this  research  was  accomplished  by  taking  various 
oxidizer/fuel  combinations  and  determining  the  theoretical  performance  parameters  in 
order  to  use  them  in  future  design  of  hybrid  rocket  motors.  The  goal  was  to  determine 
the  influence  of  equivalence  ratio  on  combustion  temperature,  specific  impulse  (Isp),  ratio 
of  specific  heats  of  combustion  products  (y),  and  molecular  mass  of  combustion  products. 
Fuel/oxidizer  combinations  analyzed  were  50%/50%  (by  mass)  lard/paraffin  and  O2  (g), 
50%/50%  lard/paraffin  and  N20  (g),  beeswax  and  02  (g),  beeswax  and  N20  (g),  lard  and 
02  (g),  lard  and  N20  (g),  beeswax  +  A1  and  02  (g),  and  beeswax  +  A1  and  N20  (g). 
Additives  to  the  solid  fuel  were  examined,  such  as  powdered  aluminum,  to  increase  the 
density  of  the  fuel  and  the  specific  impulse.  For  the  beeswax/Al  combinations,  various 
percentages  of  aluminum  were  examined. 

The  analysis  was  accomplished  using  the  Air  Force  Chemical  Equilibrium 
Specific  Impulse  Code,  also  referred  to  as  the  Isp  code  (23).  It  has  an  input  file  utility 
written  by  Curtis  Selph  which  aids  in  setting  up  the  inputs  to  the  program.  This  code 
accomplishes  the  thermodynamic  calculations  required  to  determine  a  theoretical  value 
for  the  specific  impulse  (Isp)  for  a  fucl/oxidizcr  combination.  It  takes  as  inputs  the 
propellants,  the  proportions  of  the  propellants,  and  the  combustion  chamber  pressure.  It 
can  be  specified  either  to  set  an  area  ratio  for  the  nozzle  or  to  set  a  specific  nozzle  exit 
pressure  in  which  case  the  code  will  compute  the  required  nozzle  area  ratio.  When 
calculating  the  conditions  in  the  combustion  chamber,  the  code  accounts  for  all  possible 
chemical  species  during  the  combustion  and  accounts  for  dissociation.  An  example  of 
the  mole  fractions  of  different  chemical  species  present  in  the  combustion  process  of 
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beeswax  and  gaseous  oxygen  is  shown  in  Table  1.  These  were  computed  at  a  combustion 
chamber  pressure  of  500.38  psig  and  at  an  equivalence  ratio  of  0.8.  There  are  two  modes 
for  calculating  the  flow  of  the  combustion  products  through  the  nozzle.  By  default,  the 
code  runs  performs  the  calculations  by  assuming  the  flow  through  the  nozzle  is  in  a 
shifting  equilibrium.  This  assumes  that  at  any  point  in  the  nozzle,  the  flow  is  always  in 
equilibrium  as  it  expands  through  the  nozzle  and  the  temperature  and  pressure  are 
changing.  This  tends  to  over  predict  the  Isp  computed  because  the  flow  generally  will 
travel  too  quickly  through  the  nozzle  as  the  temperature  decreases  to  be  able  to  complete 
all  the  chemical  reactions  in  order  to  be  in  chemical  equilibrium.  The  other  option  is  to 
choose  the  frozen  flow  mode  of  calculation.  This  assumes  the  chemical  composition  is 
static  as  the  gas  flow  expands  through  the  nozzle.  The  gas  can  be  designated  as  frozen  at 
the  nozzle  inlet,  the  nozzle  throat,  or  the  nozzle  exit.  For  this  research,  the  gas  flow  is 
assumed  to  be  frozen  at  the  combustion  chamber.  Although  it  is  common  to  assume 
frozen  flow  at  the  throat,  there  is  little  difference  in  assuming  frozen  flow  at  the  throat  or 
the  nozzle  inlet.  However,  setting  frozen  flow  at  the  throat  requires  the  nozzle  area  ratio 
to  be  set  in  the  input  file.  This  is  undesirable  since  the  experimental  data  that  was 
compared  with  the  theoretical  data  was  adjusted  to  not  include  any  pressure  thrust  effects, 
which  means  the  experimental  data  is  for  a  nozzle  expanded  to  sea  level  pressure  at  the 
exit.  It  is  easier  to  let  the  code  calculate  the  required  nozzle  area  ratio,  especially  when 
the  combustion  chamber  pressure  is  changed. 

The  input  file  was  setup  to  evaluate  the  Isp  with  the  nozzle  exit  pressure  set  at  14.7 
psig.  The  code  was  run  at  two  different  combustion  pressures.  The  first  combustion 
pressure  examined  was  150  psig.  This  pressure  was  chosen  because  it  was  approximately 
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Table  1:  Combustion  Chemical  Species  Mole  Fractions,  Beeswax/02,  a  =  0.8 


Species 

Mole  Fraction 

H 

0.038 

h2 

0.076 

HO 

0.078 

ho2 

0.00017 

h2o 

0.307 

h2o2 

0.000012 

0 

0.022 

o2 

0.035 

HCOOH 

0.0000023 

CHO 

0.000016 

CO 

0.297 

C02 

0.145 
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in  the  middle  of  the  range  of  the  experimental  combustion  chamber  pressures  achieved 
during  testing.  This  allowed  for  a  valid  comparison  between  the  theoretical  data  and  the 
experimental  data.  The  second  combustion  pressure  chosen  was  500.38  psig  (3.45  MPa). 
This  allows  for  comparison  to  previously  published  data  for  traditional  fuel/oxidizer 
combinations,  such  as  in  Reference  1 . 

The  Isp  code  has  a  built  in  library  of  propellants  from  which  the  user  may  choose. 
The  oxidizers  of  interest,  oxygen  and  nitrous  oxide,  already  existed  in  the  library,  but 
lard,  paraffin,  and  beeswax  were  not  in  the  library.  These  propellants  had  to  be  entered 
into  the  library  manually.  The  chemical  formula,  heat  of  formation,  and  density  and 
standard  temperature  and  pressure  were  entered  into  the  library  for  each  propellant.  Lard 
and  paraffin  were  straightforward  to  enter  since  the  properties  are  easy  to  find  (24). 
Beeswax  was  not  so  easy  to  find  exact  information  on,  so  extra  effort  was  required  to 
produce  the  information  needed.  Beeswax  is  a  complex  substance  used  in  beehives  for 
the  construction  of  honeycomb  (25).  The  worker  bees  in  the  hive  secrete  beeswax  from 
their  abdomens  and  then  combined  with  propolis  to  make  the  honeycomb.  Propolis  is  a 
substance  made  of  resins  that  bees  extract  from  the  buds  of  trees.  Beeswax  is  a  complex 
mixture  of  various  esters,  hydrocarbons,  free  acids,  alcohols,  and  other  substances  (26- 
28).  The  exact  chemical  composition  is  dependent  on  what  type  of  bee  produces  the 
beeswax  and  the  type  of  pollen  consumed  by  the  bees.  In  general,  most  beeswax  is  made 
up  primarily  of  three  main  constituents.  Cerotic  acid  (CH3(CH2)24COOH)  is  the  most 
prevalent  free  acid  and  makes  up  approximately  12%  of  beeswax.  Triacontanol 
(CH3(CH2)29CH3)  is  an  ester  that  makes  up  roughly  65%  of  beeswax.  The  third  major 
component  of  beeswax  is  myricyl  pahnitate  (CH3(CH2)i4COO(CH2)i2)CH3),  which  is 
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approximately  23%  of  beeswax.  Using  these  three  components  and  their  respective 
contributions,  a  general  chemical  formula  can  be  estimated  to  be  C46H92O.  The  data  for 
lard,  beeswax,  and  paraffin  needed  for  use  in  the  code  is  listed  in  Table  2  (24).  The  range 
of  equivalence  ratios  examined  was  0.2  -  3.  This  was  achieved  by  running  the  propellant 
combinations  at  a  series  of  oxidizer  to  fuel  ratios  that  resulted  in  the  desired  range  of 
equivalence  ratios  base  on  the  stoichiometric  oxidizer  to  fuel  ratios.  The  stoichiometric 
oxidizer  to  fuel  ratios  for  the  fuel/oxidizer  combinations  are  listed  in  Table  3. 

The  stoichiometric  ratios  were  computed  by  balancing  the  fuel  components  with 
the  exhaust  products  in  the  chemical  formulas  for  complete  combustion  and  then 
performing  an  oxygen  balance  on  each  equation. 


Table  2:  Physical  Data  on  Solid  Fuels  at  Standard  Temperature  (24) 


Fuel 

Chemical  Formula 

Heat  of  Formation  (kcaPmol) 

Density  (g/cc) 

Beeswax 

C46H92O 

-197.858 

0.961 

Lard 

C15H26O6 

-249.9 

0.649 

Paraffin 

C20H42 

-110.42 

0.9 
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Table  3:  Stoichiometric  Oxidizer  to  Fuel  Ratios  By  Mass 


Propellants 

Stoichiometric  Oxidizer  to  Fuel  Ratio 

Beeswax  +  Oxygen 

3.32 

Beeswax  +  Nitrous  Oxide 

9.13 

Lard  +  Oxygen 

1.96 

Lard  +  Nitrous  Oxide 

5.39 

50/50  Lard/Paraffin  +  Oxygen 

2.71 

50/50  Lard/Paraffin  +  Nitrous  Oxide 

7.45 

Beeswax  +  5%  A1  +  Oxygen 

3.19 

Beeswax  +  10%  A1  +  Oxygen 

3.08 

Beeswax  +  15%  A1  +  Oxygen 

2.96 

Beeswax  +  20%  A1  +  Oxygen 

2.83 

The  following  chemical  equations  for  the  combustion  were  used: 

C^6H920(s)  +  68.502(g)  — >  46C02(g)  +  46H20(g)  (4) 

C46H920(s)  +  \31  N  20(g)  -+  46C02  (g)  +  46 H20(g)  + 137 N2  (g)  (5) 

C15H2606  (s)  + 1 8.502 (g)  ->  1 5C6>2  (g)  + 13 H20(g)  (6) 

Cl5H2606 (s)  +  37 N20(g)  ->  15 C02  (g)  + 13 H20(g)  +  37 N2  (g)  (7) 

C15H2606(s)  +  1.07 C20H42(s)  +  51.33 02(g)  ->  36.53C02(g)  +  35.6 H20(g)  (8) 


Cl5//260(1(.s)+ 1 . 0  7C7  ()  /  /4  2  (.v )  + 1  °265.'V20(g)  -+36.5X'02(g)+35.6//2fTg)  +  l  °265A''2(g)  (9) 

Eqs.  4  and  5  represent  the  combustion  of  beeswax  with  oxygen  and  nitrous  oxide, 
respectively.  Similarly,  Eqs.  6  and  7  represent  the  combustion  of  lard.  Eqs.  8  and  9 
illustrate  the  combustion  of  a  mixture  of  50%  lard  and  50%  paraffin.  For  the  addition  of 
aluminum,  various  percentages  of  aluminum  by  weight  were  added  to  the  beeswax.  A 
range  of  5%  -  20%  of  aluminum  by  mass  was  chosen  for  analysis  (2). 
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III.  ANALYTICAL  RESULTS 


III.I  Results  for  =  150  psig 

The  first  set  of  analytical  data  was  examined  at  a  combustion  chamber  pressure  of 
150  psig  in  order  to  allow  for  comparison  with  the  experimental  data.  Individual  plots  of 
temperature,  specific  impulse,  ratio  of  specific  heat,  and  combustion  products  molecular 
mass  vs.  equivalence  ratio  are  given  in  Appendix  A.  Table  4  below  shows  the  maximum 
temperature  and  vacuum  specific  impulse  for  each  of  the  fuel/oxidizer  combinations  that 
were  analyzed.  Vacuum  specific  impulse  is  the  resulting  Isp  if  the  rocket  was  exhausting 
into  a  vacuum  as  opposed  to  the  atmosphere.  There  are  several  conclusions  that  can  be 
gleaned  from  Table  4.  Beeswax  with  the  addition  of  aluminum  has  the  highest  maximum 
temperature,  with  the  temperature  increasing  as  the  percentage  of  aluminum  increases. 
Beeswax,  the  lard/paraffin  mixture,  and  lard  with  oxygen  have  the  next  highest 
temperatures,  in  that  order.  It  is  worth  noting  that  adding  20%  aluminum  by  mass  only 
increases  the  maximum  temperature  by  1 18  K  over  no  aluminum.  Beeswax,  the 
lard/paraffin  mixture,  and  lard  with  nitrous  oxide  have  the  next  lowest  maximum 
temperatures,  respectively.  These  relationships  can  be  seen  in  Figure  2.  It  can  be 
determined  from  Figure  2  that  the  maximum  temperature  occurs  at  an  equivalence  ratio 
of  approximately  0.8. 

It  should  be  noted  that  the  temperatures  listed  in  Table  4  are  the  temperatures  in 
the  combustion  chamber  and  not  at  the  nozzle  exit.  The  temperature  at  the  nozzle  exit 
will  be  much  less  than  the  chamber  temperature  due  to  expansion  through  the  nozzle.  For 
example,  the  combustion  chamber  temperature  for  a  typical  equivalence  ratio  (i.e.  0.3)  in 
the  experimental  testing  would  be  approximately  1657  K  as  read  from  Figure  2. 
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Table  4:  Temperature  and  Specific  Impulse  Results,  Pc  =  150  psig 


Propellants 

Maximum  Temperature  (K) 

Maximum  Isp  (s) 

Beeswax  +  Gaseous  Oxygen 

3462 

267 

Lard  +  Gaseous  Oxygen 

3355 

250 

50/50  Lard/Paraffin  +  Gaseous 

Oxygen 

3416 

260 

Beeswax  +  Nitrous  Oxide 

3237 

237 

Lard  +  Nitrous  Oxide 

3168 

231 

50/50  Lard/Paraffin  +  Nitrous  Oxide 

3210 

234 

Beeswax  +  5%  Al  +  Gaseous  Oxygen 

3482 

266 

Beeswax  +  10%  Al  +  Gaseous 

Oxygen 

3518 

266 

Beeswax  +  15%  Al  +  Gaseous 

Oxygen 

3549 

268 

Beeswax  +  20%  Al  +  Gaseous 

Oxygen 

3580 

269 

Paraffin  +  Gaseous  Oxygen 

3447 

223 

Temperature  vs  Equivalence  Ratio 
Combustion  Chamber  Pressure  =  150  psig 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8  2  2.2  2.4  2.6  2.8  3 

Equivalence  Ratio 


Beeswax  +  Oxygen 
^“•Beeswax  +  Nitrous  Oxide 
Lard  +  Oxygen 
Lard  +  Nitrous  Oxide 
^“50/50  Lard/Paraffin  +  Oxygen 
^“50/50  Lard/Paraffin  +  Nitrous  Oxide 
Beeswax  and  5%  Al  +  Oxygen 
^“Beeswax  and  10%  Al  +  Oxygen 
Beeswax  and  15%  Al  +  Oxygen 
Beeswax  and  20%  Al  +  Oxygen 


Figure  2:  Comparison  of  Temperature  vs.  Equivalence  Ratios,  Pc  =  150  psig 
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At  this  same  equivalence  ratio,  the  ratio  of  specific  heats  is  approximately  1.3  from 
Appendix  A.  Based  on  the  expansion  ratio  of  the  nozzle,  the  exit  plane  Mach  number  is 
approximately  2.3.  Using  these  parameters  and  assuming  frozen  flow  in  the  nozzle,  the 
ratio  of  combustion  chamber  temperature  to  exit  plane  temperature  can  be  calculated 
fromEq.  10: 


Zo 

Te 


(10) 


where  T0  is  the  combustion  chamber  temperature,  Te  is  the  exit  temperature,  Me  is  the  exit 
plane  Mach  number,  and  y  is  the  ratio  of  specific  heats  (29).  In  this  example,  the 
resultant  ratio  of  chamber  temperature  to  exit  temperature  is  1.8.  This  yields  an  exit 
plane  temperature  of  approximately  920  K. 

Figure  3  shows  a  similar  comparison  between  the  fucl/oxidizer  combinations  for 
vacuum  specific  impulse.  It  is  interesting  to  note  the  behavior  when  aluminum  is  added 
to  beeswax  and  burned  with  oxygen.  The  Isp  is  virtually  the  same  with  no  aluminum 
added  as  with  any  level  of  aluminum  addition  considered.  Previous  analysis  has  shown 
that  addition  of  aluminum  to  traditional  hybrid  solid  fuels  actually  decreases  the  specific 
impulse  because  the  molecular  weight  of  the  exhaust  products  is  increased,  which  negates 
the  increase  in  the  temperature  (2). 

When  determining  the  specific  impulse  of  an  ideally  expanded  rocket  motor,  the 
main  driver  is  the  exhaust  velocity.  Eq.  1 1  shows  the  dependence  of  the  exit  velocity  on 
the  chamber  temperature  and  the  molecular  mass  of  the  combustion  products 
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Vacuum  Specific  Impulse  vs  Equivalence  Ratio 
(Frozen  Flow,  Pc  =  150  psig) 


^—Beeswax  +  Oxygen 
““Beeswax  +  Nitrous  Oxide 
Lard  +  Oxygen 
Lard  +  Nitrous  Oxide 
““50/50  Lard/Paraffin  +  Oxygen 
““50/50  Lard/Paraffin  +  Nitrous  Oxide 
““Beeswax  and  5%  Al  +  Oxygen 
^“Beeswax  and  10%  Al  +  Oxygen 
^“Beeswax  and  15%  Al  +  Oxygen 
^“Beeswax  and  20%  Al  +  Oxygen 


Equivalence  Ratio 


Figure  3:  Comparison  of  Vacuum  Specific  Impulse  vs.  Equivalence  Ratio,  Pc  =  150  psig 
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where  Ve  is  the  nozzle  exit  velocity,  y  is  the  ratio  of  specific  heats,  To  is  the  combustion 
chamber  temperature,  Mis  the  combustion  products  average  molecular  mass,  Ru  is  the 
universal  gas  constant,  Pe  is  the  nozzle  exit  plane  pressure,  and  Po  is  the  combustion 
chamber  pressure.  For  our  case,  the  maximum  Isp  occurs  at  an  equivalence  ratio  of 
approximately  0.6.  This  is  roughly  the  optimal  equivalence  ratio  to  operate  a  rocket 
motor  (1).  It  is  interesting  to  note  that  as  more  aluminum  is  added,  the  equivalence  ratio 
where  the  maximum  Isp  occurs  tends  to  shift  to  the  left.  The  reason  the  equivalence  ratio 
where  the  maximum  Isp  occurs  is  shifted  to  the  left  as  compared  to  where  the  maximum 
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chamber  temperature  occurs  is  due  to  Eq.  1 1 .  Eq.  1 1  is  a  function  of  chamber 
temperature,  molecular  mass  of  the  combustion  products,  and  the  ratio  of  specific  heats. 
Figures  4  and  5  show  the  behavior  of  the  average  molecular  mass  of  the  combustion 
products  and  the  ratio  of  specific  heats  as  a  function  of  the  equivalence  ratio  for  beeswax 
and  oxygen,  respectively.  The  molecular  mass  clearly  increases  as  the  equivalence  ratio 
increases.  The  ratio  of  specific  heats  decreases  as  the  equivalence  ratio  increases. 
However,  the  change  in  the  absolute  value  is  not  significant  enough  to  overcome  other 
terms  in  Eq.  11.  It  is  the  combination  of  these  factors  that  shifts  the  equivalence  ratio  to 
the  left  for  the  maximum  Isp.  The  best  way  to  illustrate  this  is  to  plot  the  ratio  of  the 
combustion  temperature  to  the  molecular  mass  of  the  combustion  products.  This  is  the 
dominant  term  in  Eq.  1 1  for  determining  where  the  maximum  Isp  occurs. 


Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  150  psig) 


Figure  4:  Molecular  Mass  vs.  Equivalence  Ratio,  Pc  =  150  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  150  psig) 


Figure  5:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Pc  =  150  psig 


As  can  be  seen  from  Figure  6,  the  maximum  of  the  ratio  of  temperature  to  molecular 
mass  occurs  at  an  equivalence  ratio  of  approximately  0.6.  This  is  exactly  where  the  Isp 
curve  is  at  its  maximum. 

From  Figure  7,  we  can  see  the  behavior  of  the  ratio  of  specific  heats  vs 
equivalence  ratio  for  various  propellant  combinations.  The  combinations  using  nitrous 
oxide  have  the  highest  values.  Several  trends  are  noticeable  in  the  figure.  As  the 
percentage  of  aluminum  increases,  the  ratio  of  specific  heats  decreases.  The  curve  of  the 
lard/paraffin  mixture  with  oxygen  is  almost  identical  to  the  beeswax  plus  5%  aluminum 
with  oxygen  curve.  The  curve  for  lard  plus  oxygen  is  also  very  close  to  the  curve  for 
beeswax  plus  10%  aluminum  with  oxygen. 
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Ratio  of  Specific  Heats 


Ratio  of  Temperature  to  Combustion  Product  Molecular 
Weight  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  150  psig) 


Figure  6:  Ratio  of  Temperature  to  Molecular  Mass  vs.  Equivalence  Ratio,  Pc  =  150  psig 


Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
(Pc  =  150  psig) 


- Beeswax  +  Oxygen 

- Beeswax  +  Nitrous  Oxide 

- Lard  +  Oxygen 

Lard  +  Nitrous  Oxide 

- 50/50  Lard/Paraffin  +  Oxygen 

- 50/50  Lard/Paraffin  +  Nitrous  Oxide 

- Beeswax  and  5%  Al  +  Oxygen 

- Beeswax  and  10%  Al  +  Oxygen 

- Beeswax  and  15%  Al  +  Oxygen 

Beeswax  and  20%  Al  +  Oxygen 


Figure  7:  Comparison  of  Ratio  of  Specific  Heats,  Pc  =  150  psig 
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In  Figure  8,  the  average  molecular  mass  of  the  combustion  products  are  shown  versus 
equivalence  ratio  for  various  propellant  combinations.  As  expected,  the  molecular  mass 
increases  as  the  addition  of  aluminum  to  beeswax  increases  from  zero  to  20%.  The  three 
curves  for  nitrous  oxide  as  the  oxidizer  behave  somewhat  differently  than  those  for 
oxygen.  These  curves  seem  to  be  asymptotically  approaching  a  value  between  29  and  30 
g/mol  as  the  equivalence  ratio  increases.  It  is  unclear  from  the  range  of  equivalence 
ratios  studied  here  if  the  molecular  mass  of  the  fuel/oxidizer  combinations  using  oxygen 
will  asymptotically  approach  some  value  as  the  equivalence  ratio  continues  to  increase. 


Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
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- Beeswax  +  Oxygen 

- Beeswax  +  Nitrous  Oxide 

- Lard  +  Oxygen 

- Lard  +  Nitrous  Oxide 

- 50/50  Lard/Paraffin  +  Oxygen 

- 50/50  Lard/Paraffin  +  Nitrous  Oxide 

- Beeswax  and  5%  Al  +  Oxygen 
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- Beeswax  and  20%  Al  +  Oxygen 


Equivalence  Ratio 


Figure  8:  Comparison  of  Molecular  Mass,  Pc  =  150  psig 
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III.II  Results  for  =  500.38  psig 

A  second  set  of  analytical  data  was  calculated  for  a  combustion  chamber  pressure 
of  500.38  psig  (3.45  MPa)  to  allow  for  comparison  to  previously  published  data  for  more 
traditional  hybrid  rocket  motors.  In  addition  to  the  change  in  combustion  chamber 
pressure,  the  nozzle  exit  pressure  was  changed  to  2.0  psig  to  match  the  conditions  that  the 
data  in  Reference  1  were  calculated  for.  Appendix  B  contains  plots  for  temperature, 
vacuum  specific  impulse,  ratio  of  specific  heats,  and  molecular  mass  of  the  combustion 
products.  The  traditional  hybrid  propellant  combinations  used  for  comparison  were 
HTPB/hydrogen  peroxide,  HTPB/LOX,  and  HTPB/nitrogen  tetroxide.  These  three 
oxidizers  are  good  oxidizers  but  have  significant  safety  and  handling  issues  (2).  The 
hydrogen  peroxide  considered  is  90  %  hydrogen  peroxide  and  10%  water.  Hydrogen 
peroxide  is  toxic  and  a  fire  hazard.  LOX  is  cryogenic,  not  easily  storable,  difficult  to 
handle,  and  an  explosive  hazard.  Nitrogen  tetroxide  is  the  most  common  storable 
oxidizer  in  the  U.S.,  but  has  toxic  fumes. 

A  limited  validation  study  was  performed  to  evaluate  the  accuracy  of  the  Isp  code. 
Three  propellant  combinations  were  examined:  kerosene/LOX,  liquid  hydrogen/LOX, 
and  UDMH/LOX.  The  maximum  temperature  and  sea  level  specific  impulse  were  found 
at  a  combustion  chamber  pressure  of  1000  psig  from  Reference  30.  The  Isp  code  was 
used  to  calculate  the  maximum  temperature  and  sea  level  specific  impulse  at  the  same 
combustion  chamber  pressure.  A  summary  of  the  results  can  be  seen  in  Table  5.  The 
largest  difference  (3.8%)  was  in  the  maximum  temperature  of  liquid  hydrogen/LOX. 
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Table  5:  Known  vs.  Calculated  Maximum  Temperatures  and  Sea  Level  Isp 


Propellant  Combination 

Kerosene/LOX 

Liquid  Hydrogen/LOX 

UMDH//LOX 

Optimal  O/F  Ratio 

2.29 

5.0 

1.38 

Known  Max.  Temperature 
(K) 

3540 

3425 

3460 

Isp  Code  Max. 

Temperature  (K) 

3577 

3295 

3460 

%  Difference  in 
Temperature 

1.0 

3.8 

0.4 

Known  Max.  Isp  (sec) 

289 

381 

297 

Isp  Code  Max.  Isp  (sec) 

285 

378 

299 

%  Difference  in  Isp 

1.4 

0.8 

0.7 
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The  maximum  temperatures  and  vacuum  specific  impulses  for  the  analyzed 
fuel/oxidizer  combinations  are  shown  in  Table  6.  The  same  trends  are  evident  in  the  data 
at  500.38  psig  as  in  the  data  at  150  psig.  The  average  temperature  increase  due  to 
increasing  the  chamber  pressure  from  150  psig  to  500.38  psig  was  168.9  K.  The 
corresponding  average  increase  in  the  vacuum  specific  impulse  was  52.5  s.  These  are 
relatively  small  increases  when  compared  to  the  large  increase  in  the  combustion 
chamber  pressure.  The  comparison  of  the  temperatures  in  Table  6  is  shown  in  Figure  9. 
The  propellant  combinations  using  oxygen  as  the  oxidizer  resulted  in  higher  temperatures 
than  when  using  nitrous  oxide.  Regardless  of  the  oxidizer,  beeswax  produced  the  highest 
temperature  and  lard  produced  the  lowest  except  for  HTPB/hydrogen  peroxide  and 
HTPB/nitrogen  tetroxide.  As  the  amount  of  aluminum  added  to  the  beeswax  increases, 
the  temperature  increases.  Similarly  to  the  data  at  a  pressure  of  150  psig,  the  maximum 
temperature  at  500.38  psig  occurs  at  an  equivalence  ratio  of  approximately  0.8.  Also 
shown  in  Table  6  are  other  propellant  combinations.  Maximum  temperatures  for  liquid 
oxygen  combined  with  kerosene  and  liquid  hydrogen  are  shown  as  taken  from  reference 
1.  Also  shown  is  the  maximum  temperature  of  paraffin  with  gaseous  oxygen,  beeswax 
with  liquid  oxygen,  and  UDMH  with  liquid  oxygen  as  computed  using  the  Air  Force 
Specific  Impulse  Code.  These  temperatures  are  all  in  the  range  of  3500  K  to  3630  K,  and 
are  very  similar  to  the  temperatures  given  for  the  non-conventional  propellants  analyzed 
in  this  study.  These  results  support  the  validity  of  the  maximum  combustion  chamber 
temperatures  presented  in  Table  5  for  the  propellant  combinations  considered  here. 

The  results  for  the  vacuum  specific  impulse  behave  similarly  to  the  temperature, 
except  the  maximum  occurs  at  an  equivalence  ratio  of  0.6.  The  propellant  combinations 
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Table  6:  Temperature  and  Specific  Impulse  Results,  Pc  =  500.38  psig 


Propellants 

Maximum 
Temperature  (K) 

Maximum  Isp  (s) 

Beeswax  +  Gaseous  Oxygen 

3650 

327 

Lard  +  Gaseous  Oxygen 

3524 

308 

50/50  Lard/Paraffm  +  Gaseous 

Oxygen 

3600 

303 

Beeswax  +  Nitrous  Oxide 

3365 

286 

Lard  +  Nitrous  Oxide 

3288 

280 

50/50  Lard/Paraffin  +  Nitrous  Oxide 

3328 

271 

Beeswax  +  5%  Al  +  Gaseous  Oxygen 

3670 

312 

Beeswax  +  10%  Al  +  Gaseous 

Oxygen 

3711 

328 

Beeswax  +  15%  Al  +  Gaseous 

Oxygen 

3747 

328 

Beeswax  +  20%  Al  +  Gaseous 

Oxygen 

3783 

330 

Paraffin  +  Gaseous  Oxygen 

3630 

307 

Beeswax  +  Liquid  Oxygen 

3606 

326 

Traditional  Propellants 

HTPB  +  Hydrogen  Peroxide 

2802 

298 

HTPB  +  Liquid  Oxygen 

3669 

318 

HTPB  +  Nitrogen  Tetroxide 

3450 

297 

Temperature  vs  Equivalence  Ratio 
(Pc  =  500.38  psig) 


Equivalence  Ratio 


Beeswax  +  Oxygen 
““  Beeswax  +  Nitrous  Oxide 
^^“Lard  +  Oxygen 
^—Lard  +  Nitrous  Oxide 
““50/50  Lard/Paraffin  +  Oxygen 
^“50/50  Lard/Paraffin  +  Nitrous  Oxide 
““  Beeswax  and  5%  Al  +  Oxygen 
^“Beeswax  and  10%  Al  +  Oxygen 
Beeswax  and  1 5%  Al  +  Oxygen 
-  Beeswax  and  20%  Al  +  Oxygen 


Figure  9:  Temperature  vs.  Equivalence  Ratio,  Pc  =  500.38  psig 
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using  oxygen  had  higher  vacuum  specific  impulse  values  than  when  using  nitrous  oxide. 
Beeswax  had  the  highest  vacuum  specific  impulse  followed  by  the  lard/paraffin  mixture, 
and  then  lard.  When  5%  aluminum  was  added  to  the  beeswax,  it  reduced  the  specific 
impulse  by  15  s.  When  the  amount  of  aluminum  is  increased  to  10%  and  15%,  the 
specific  impulse  rises  to  1  s  above  the  value  achieved  with  no  aluminum.  At  20% 
aluminum,  the  specific  impulse  rises  to  3  s  about  the  value  achieved  with  no  aluminum. 
These  trends  can  be  seen  in  Figure  10. 

While  the  contribution  of  aluminum  shows  the  potential  for  slight  improvements 
in  the  vacuum  specific  impulse,  the  negative  effects  of  using  aluminum  would  likely 
negate  this  advantage.  When  using  an  aluminized  propellant,  such  as  in  solid  rocket 
motors,  there  is  significant  nozzle  erosion  due  to  impingement  of  aluminum  oxide 
particles  with  the  nozzle  walls.  If  severe  enough,  this  erosion  can  cause  failure  of  the 


Vacuum  Specific  Impulse  vs  Equivalence  Ratio 
(Pc  =  500.38  psig) 


^“Beeswax  +  Oxygen 
•^Beeswax  +  Nitrous  Oxide 
Lard  +  Oxygen 
Lard  +  Nitrous  Oxide 
““50/50  Lard/Paraffin  +  Oxygen 
““50/50  Lard/Paraffin  +  Nitrous  Oxide 
^“Beeswax  and  5%  Al  +  Oxygen 
^“Beeswax  and  10%  Al  +  Oxygen 
^“Beeswax  and  15%  Al  +  Oxygen 
““Beeswax  and  20%  Al  +  Oxygen 


Equivalence  Ratio 


Figure  10:  Vacuum  Specific  Impulse  vs.  Equivalence  Ratio,  Pc  =  500.38  psig 
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nozzle.  Another  effect  is  that  using  an  aluminized  propellant  usually  results  in  two-phase 
flow  in  the  nozzle.  As  the  gas  flow  expands  through  the  nozzle,  the  temperature  of  the 
gas  can  decrease  considerably  depending  on  the  expansion  ratio  of  the  nozzle,  which 
could  results  in  liquid  droplets  condensing  out  of  the  gas  flow.  The  presence  of  liquid 
droplets  in  the  exhaust  flow  reduces  the  specific  impulse  of  the  rocket.  This  is  due  to  the 
energy  available  to  accelerate  the  gas  flow  being  absorbed  by  the  heavier  liquid  droplets. 
These  effects  are  not  reproduced  by  the  Isp  code,  but  should  be  considered  when  adding 
aluminum  to  a  propellant. 

The  HTPB/LOX  combination  has  a  higher  maximum  temperature  than  all  of  the 
propellant  combinations  analyzed  except  when  aluminum  is  added  to  beeswax.  The 
difference  between  the  maximum  temperature  for  HTPB/LOX  and  beeswax/oxygen  is 
only  19  K.  However,  the  maximum  specific  impulse  for  beeswax/oxygen  is  9  s  higher 
than  for  HTPB/LOX.  Beeswax/oxygen  produces  higher  chamber  temperatures  and 
specific  impulses  than  either  HTPB/hydrogen  peroxide  or  HTPB/nitrogen  tetroxide. 

When  comparing  the  propellant  combinations  using  nitrous  oxide  as  the  oxidizer,  all 
three  have  a  lower  maximum  temperature  than  HTPB/LOX  and  HTPB/nitrogen  tetroxide. 
All  three  combinations  have  a  higher  temperature  than  HTPB/hydrogen  peroxide.  When 
the  maximum  vacuum  specific  impulse  is  compared,  all  the  propellant  combinations 
using  nitrous  oxide  produced  lower  specific  impulses  than  the  three  traditional  propellant 
combinations  using  HTPB. 

These  results  indicate  that  beeswax/oxygen  should  produce  the  best  performing 
hybrid  rocket  motor  when  compared  to  the  traditional  hybrid  rocket  propellants.  While 
the  performance  between  beeswax/oxygen  is  very  close  to  HTPB/LOX,  using  gaseous 
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oxygen  is  much  more  convenient  than  using  LOX.  There  are  several  factors  to  consider 
when  evaluating  what  propellants  to  use:  performance,  toxicity,  ease  of  handling,  and 
public  acceptability.  Given  the  small  performance  differences  between  the  propellants 
considered  here,  beeswax  and  gaseous  oxygen  compare  well  as  potentially  new  hybrid 
rocket  motor  propellants. 


32 


IV.  DESCRIPTION  OF  EXPERIMENTAL  PROGRAM 


IV.I  Description  of  Test  Facility 

The  hybrid  rocket  test  facility  is  located  on  the  University  of  Tennessee 
Agricultural  Research  Land  site  just  south  of  the  University  of  Tennessee  campus  in 
Knoxville,  Tennessee.  The  test  facility  consists  of  a  test  stand,  a  storage  building,  and  a 
safety  wall.  Figure  1 1  illustrates  the  setup  of  the  test  site.  The  safety  wall  was 
constructed  to  act  as  a  barrier  to  protect  individuals  in  case  of  a  mishap.  The  wall  is 
constructed  of  pressure  treated  wood  and  filled  with  sand.  Its  dimensions  are  ten  ft.  long 
by  two  and  a  half  ft.  thick  by  six  ft.  high.  The  interior  of  the  wall  contains  approximately 
six  tons  of  sand.  The  storage  building  is  used  to  securely  store  the  test  stand,  power 
generator,  oxidizer  and  propane  tanks,  and  other  test  equipment.  Due  to  the  remote 
nature  of  the  test  site,  a  generator  is  used  to  provide  electricity  for  the  instrumentation, 
electric  solenoids,  and  computers. 

The  test  stand  consists  of  a  support  cage  containing  the  rocket  motor  mounted  on 
a  transportable  steel  table  as  illustrated  in  Figure  12.  The  solid  fuel  is  formed  in  a  steel 
cylinder  which  has  a  3  in.  outer  diameter,  is  0.25  in.  thick,  and  is  10  in.  long.  The 
combustion  port  diameter  of  all  motors  tested  was  0.5  in.  The  cylinder  is  held  in  place  by 
six,  3/8  in.  diameter  bolts.  At  the  inlet  end,  there  is  a  steel  injector  end  cap  which  houses 
the  spark  plug  for  ignition.  The  oxidizer  and  propane  lines  both  feed  into  the  steel  inlet 
end  cap.  When  the  test  is  initiated,  solenoid  valves  on  the  propane  and  oxidizer  lines  are 
both  opened,  and  the  spark  plug  is  fired.  A  flow  diagram  of  the  test  setup  can  be  seen  in 
Figure  13. 
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Figure  11:  Test  Site  Setup 
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Figure  12:  Test  Stand 
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Orifice  Shutoff  Valve 


X*  X3  ®  O 

Regulator  Solenoid  Valve  Hash  Arrestor  Pressure  Transducer 


Figure  13:  Flow  Diagram  of  Flybrid  Rocket  Motor 


35 


The  propane/oxidizer  mixture  is  ignited  by  the  spark  plug  and  provides  the  initial  heat  of 
vaporization  to  the  fuel.  The  propane  solenoid  valve  and  spark  plug  are  both  controlled 
by  a  single  timer,  and  they  are  typically  on  for  0.5  s.  The  oxidizer  solenoid  is  controlled 
by  a  separate  timer,  which  typically  runs  for  5  to  7  s.  At  the  end  of  the  test,  the  timer 
automatically  closes  the  oxidizer  valve.  The  propane  is  supplied  at  80  psig,  and  the 
oxidizer  is  supplied  in  a  range  of  200-500  psig  by  a  regulator.  The  nozzle  end  cap  houses 
a  converging-diverging  nozzle  made  of  graphite.  As  indicated  in  Figure  3,  a  strain  gage 
is  attached  to  the  aft  cantilever  support  beam  in  order  to  measure  thrust  during  the  testing. 
The  strain  gage  is  calibrated  before  every  test  to  determine  the  relationship  between  the 
voltage  output  of  the  strain  gage  and  the  thrust  produced  by  the  rocket  motor.  Calibrating 
the  strain  gage  every  test  eliminated  the  need  to  determine  the  effect  of  enviromnental 
conditions  (i.e.  ambient  temperature)  on  the  thrust  measurements.  The  calibration  was 
accomplished  by  placing  weights  in  a  bucket  which  loaded  the  strain  gage  and  then 
recording  the  corresponding  voltages  (see  Figure  14  for  an  illustration  of  the  calibration 
setup).  The  range  of  weights  used  in  the  calibration  covered  the  range  of  thrust  values 
measured  during  the  testing. 


Figure  14:  Calibration  Setup 
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Figure  15  shows  an  example  of  a  test  calibration.  The  curvefit  from  the  calibration  data 
was  used  to  convert  the  voltages  from  the  strain  gage  to  forces.  Pressure  in  the 
combustion  chamber  is  measured  using  the  pressure  transducer  positioned  downstream  of 
the  orifice  plate,  which  is  just  upstream  of  the  oxygen  inlet  into  the  combustion  chamber. 
A  restrictor  plate  was  added  after  problems  arose  with  combustion  occurring  on  the  inlet 
surface  of  the  fuel  grain,  which  led  to  burning  along  the  outside  of  the  fuel  grain.  The 
pressure  transducers  used  were  Omega  PX603s  which  have  an  accuracy  of  ±  0.4%  of  the 
full  span,  which  was  1000  psig.  Prior  experimentation  has  shown  that  the  pressure 
measured  downstream  of  the  orifice  plate  is  identical  to  the  pressure  measured  at  the 
nozzle  end  of  the  combustion  chamber  (30).  Unfortunately,  due  to  budget  restrictions, 
the  pressure  difference  across  the  orifice  plate  to  compute  the  oxidizer  mass  flow  rate  had 
to  be  determined  using  the  same  pressure  transducers. 


Calibration  Curve 
Test  169 


Figure  15:  Example  of  Calibration  Curve 
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The  nozzle  is  a  converging-diverging  type  designed  to  have  an  exit  pressure  of 
14.7  psig  for  a  combustion  chamber  pressure  of  approximately  150  psig.  For  ease  of 
fabrication,  the  nozzle  uses  a  simple  conical  design.  It  is  constructed  out  of  graphite  to 
prevent  degradation  due  to  the  high  temperature.  The  nozzle  is  2.0  in.  long  with  the 
converging  section  being  0.53  in.  long  and  a  diverging  section  that  is  1.47  in.  long.  The 
inlet,  throat,  and  exit  diameters  are  1.4  in.,  0.35  in.,  and  0.53  in.,  respectively.  The  ratio 
of  the  exit  area  to  the  throat  area  is  2.3 12,  which  results  in  an  exit  Mach  number  of 
approximately  2.3,  assuming  a  value  of  1.3  for  the  ratio  of  specific  heats.  The  goal  of  the 
nozzle  design  is  to  reduce  the  effect  of  pressure  thrust  on  the  thrust  data. 

IV.II  Mass  Flow  Rate  Uncertainty 

Table  7  shows  the  uncertainty  analysis  of  the  oxidizer  mass  flow  rate  for  test  176 
as  an  example.  The  absolute  value  of  the  range  of  uncertainties  between  the  two  worst 
case  scenarios  in  the  measurement  of  the  mass  flow  rate  is  18  -  28%.  Due  to  this 
significant  uncertainty,  the  equivalence  ratios  computed  in  this  study  also  have  a  large 
amount  of  uncertainty.  Due  to  the  large  uncertainty  in  the  oxidizer  mass  flow  rate 
measurement,  an  attempt  was  made  to  evaluate  the  uncertainty  more  accurately.  A  test 
was  run  with  only  oxidizer  flowing  through  the  system  at  an  average  upstream  pressure 
of  106  psig  for  approximately  44.5  s.  The  change  in  mass  of  the  oxidizer  tank  during  the 
test  was  1 .9  kg  ±  0. 1  kg  resulting  in  an  average  oxidizer  mass  flow  rate  of  0.043  kg/s  ± 
0.002  kg/s.  The  oxidizer  mass  flow  rate  calculated  from  the  pressure  transducer  data  was 
0.038  kg/s,  which  resulted  in  a  difference  of  1 1%  ±  4.7%.  This  indicated  that  the  actual 
uncertainty  in  the  oxidizer  mass  flow  rate  is  less  than  the  theoretical  uncertainty  due  to 
the  pressure  transducers.  The  other  source  of  uncertainty  is  the  possible  spallation  of 
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solid  particles  or  liquid  droplets  of  fuel  exiting  the  nozzle  without  being  burned.  The 
ejection  of  this  unburned  fuel  increases  the  calculated  regression  rate  without 
contributing  significantly  to  the  thrust  produced.  Table  8  shows  the  uncertainties  in  the 
instrumentation  used  in  the  experimentation.  All  of  the  instrumentation  was  wired  into  a 
Minilab  1008  A/D  converter  and  connected  to  a  laptop  via  USB  cable.  The  Minilab  1008 
has  8  channels  of  12-bit  analog  input;  however  only  three  channels  were  used:  two  for 
the  pressure  transducers  and  one  for  the  strain  gage.  The  data  was  collected  by  the  laptop 
computer  using  the  Hewlett  Packard  Visual  Engineering  Environment  (HP  VEE) 
software.  Mass  measurements  of  the  fuel  grains  were  taken  on  a  digital  scale  with  1 
gram  resolution  before  and  after  each  test. 


Table  7:  Oxidizer  Mass  Flow  Rate  Uncertainty  Estimate  for  Test  176 


Upstream 

Pressure 

(psig) 

Downstream 

Pressure 

(psig) 

Pressure 

Difference 

(psig) 

Oxidizer 
Mass  Flow 
Rate  (kg/s) 

% 

Difference 

Low  Case 

219 

210 

9 

0.0237 

28 

Experimental 

Value 

223 

206 

16 

0.0331 

N/A 

High  Case 

227 

202 

25 

0.0392 

18 

Table  8:  Uncertainties  in  Instrumentation 


Instrument 

Accuracy 

Scale 

1  gram  resolution 

Strain  Indicator 

0.1%  (full  scale) 

Timers 

1%  (set  scale) 

Pressure  Transducers 

0.4%  (full  scale) 
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IV.III  Testing  Procedure 

The  beeswax  fuel  grains  were  produced  by  melting  down  bulk  beeswax  bricks 
and  pouring  the  beeswax  into  the  steel  cylinders.  The  standard  yellow  beeswax  was 
purchased  from  the  Southeast  Texas  Honey  Company.  The  cylinders  were  placed  in 
wooden  stands  which  hold  the  cylinders  vertically  and  hold  a  polished  steel  rod  firmly  in 
the  center  of  the  cylinder.  The  steel  rod,  which  is  0.5  in.  diameter,  is  for  forming  the 
combustion  port  in  the  grain.  The  wooden  base  is  constructed  to  produce  a  space  at  the 
end  of  the  grain  1.875  in.  long  to  act  as  an  aft  mixing  chamber.  The  purpose  of  the  aft 
mixing  chamber  is  to  allow  the  gasified  propellants  to  complete  the  combustion  process 
as  much  as  possible  before  exiting  the  motor.  If  the  nozzle  entrance  were  immediately  at 
the  end  of  the  fuel  grain,  it  is  unlikely  that  the  combustion  process  would  be  complete. 
The  common  range  of  length-to-diameter  (L/D)  ratios  for  the  aft  mixing  chamber  is  0.5  - 
1.0  (1).  In  this  research,  the  aft  mixing  chamber  had  a  L/D  of  0.75.  The  liquid  beeswax 
was  poured  in  the  cylinder  to  the  top,  and  a  metal  cap  was  placed  on  top  to  hold  the  steel 
rod  in  the  center  while  the  beeswax  cooled.  The  beeswax  was  allowed  to  cool  and 
solidify  at  room  temperature.  The  cylinder  was  removed  from  the  wooden  stand,  and  the 
steel  rod  hammered  out.  The  fuel  grain  was  then  inspected  for  cracks,  and  if  cracks  are 
present,  the  grain  was  melted  down,  and  the  fabrication  process  was  repeated. 

The  testing  procedure  was  established  to  ensure  safety  of  the  experimentation  and 
to  eliminate  mistakes  that  could  result  in  a  test  failure.  A  test  supervisor,  range  safety 
supervisor,  and  calibration  supervisor  are  designated  prior  to  the  test.  The  test  supervisor 
has  overall  responsibility  for  the  safety  of  the  testing.  The  calibration  supervisor  is 
responsible  for  the  calibration  of  the  strain  gage.  The  range  safety  supervisor  is 
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responsible  for  certifying  that  the  area  is  clear  before  the  test  is  initiated.  The  test  stand  is 
rolled  out  of  the  storage  building,  and  the  wheels  are  blocked  to  prevent  the  table  from 
moving  during  a  test.  Once  the  test  stand  is  secured,  the  generator  is  started,  and  all 
electrical  systems  (instrumentation,  solenoid  valves,  computers,  etc.)  are  setup.  A  fire 
extinguisher  is  placed  on  the  personnel  side  of  the  safety  barrier  in  case  of  fire  during  a 
test.  The  fuel  grain  is  weighed,  and  the  weight  is  recorded.  The  inlet  section  is  then 
examined  for  any  traces  of  fuel  that  could  possibly  cause  a  fire  during  a  test.  If  any  fuel 
is  present,  the  section  is  disassembled  and  cleaned  before  the  test  can  be  run.  Then  the 
propane  and  oxidizer  lines  are  inspected  for  leaks.  The  flow  restrictor  plate  is  then 
installed  in  the  inlet  section.  The  restrictor  plate  has  a  hole  the  same  diameter  as  the 
combustion  port  (0.5  in.).  The  purpose  of  the  restrictor  plate  is  to  confine  the  oxidizer 
flow  to  the  combustion  port.  The  fuel  grain  is  inserted  into  the  test  stand,  and  the  bolts 
are  tightened  evenly  to  prevent  leakage  during  the  test.  At  this  point,  all  personnel  except 
for  the  test  supervisor  and  calibration  supervisor  are  evacuated  to  behind  the  safety  wall. 
The  test  supervisor  pressurizes  the  oxidizer  and  fuel  lines  by  opening  the  supply  valves. 
The  calibration  supervisor  then  calibrates  the  strain  gage.  Once  the  calibration  is 
completed,  the  test  supervisor  ensures  the  oxidizer  regulator  is  set  to  the  desired  pressure 
and  verifies  that  the  propane  and  oxidizer  timers  are  set  properly.  The  calibration 
supervisor  then  evacuates  to  the  other  side  of  the  safety  wall.  The  spark  plug  is  then 
connected,  and  the  video  recording  systems  is  started  by  the  test  supervisor.  The  test 
supervisor  moves  behind  the  safety  wall,  and  the  data  recording  system  is  activated.  The 
range  safety  supervisor  then  informs  the  test  supervisor  when  the  area  is  clear  for  the  test. 
The  firing  switch  is  plugged  in,  the  firing  switch  is  switched  from  safe  to  fire  mode,  and 
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then  the  fire  button  is  pressed  by  the  test  supervisor  to  initiate  the  test.  When  the  test 
terminates,  the  firing  switch  is  turned  back  to  safe  and  is  unplugged.  The  firing  switch  is 
only  plugged  in  for  the  duration  of  the  test  to  prevent  inadvertent  opening  of  the  solenoid 
valves.  The  video  recording  system  and  data  recording  system  are  both  stopped.  The  test 
supervisor  approaches  the  test  stand  and  disconnects  the  spark  plug.  At  this  point,  the 
test  stand  is  safe  to  approach  by  personnel.  After  the  fuel  grain  has  had  at  least  5  minutes 
to  cool  down,  it  is  removed  from  the  test  stand  and  weighed,  and  the  weight  is  recorded. 
IV.IV  Data  Analysis 

The  regression  rate  (or  burn  rate)  of  the  solid  fuel  is  defined  as  the  rate  at  which 
the  exposed  fuel  surface  is  burned  during  combustion.  It  is  measured  in  a  direction 
normal  to  the  fuel  surface  (1).  A  time-averaged  value  for  the  regression  rate  in  a  hybrid 
rocket  motor  can  be  determined  by  determining  the  change  in  the  fuel  grain  thickness  and 
divided  by  the  time  duration  of  the  test  as  shown  in  Eq.  12: 

t  —t  f 

r  =  ^ - ^  (12) 

h 

where  r  is  the  regression  rate,  tgl  is  the  initial  fuel  grain  thickness,  tg/  is  the  final  fuel 
grain  thickness,  and  tb  is  the  burn  time.  The  initial  grain  thickness  is  simply  measured 
before  a  test.  The  average  final  grain  thickness  is  determined  by  determining  the  mass  of 
the  burned  fuel  and,  using  the  known  density  of  the  fuel,  calculating  the  average 
thickness  based  on  a  change  in  volume  of  the  fuel  grain  (16). 

The  mass  flow  rate  of  the  oxidizer  was  approximated  using  a  sonic  orifice  plate 
and  two  pressure  transducers  as  shown  in  Figure  16.  The  Omega  pressure  transducers 
have  a  range  of  0-1000  psig  and  an  error  of  0.4%.  The  pressure  differential  across  the 
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orifice  is  found  by  measuring  the  pressure  just  upstream  and  immediately  downstream  of 
the  orifice  plate.  A  MATLAB  code  was  used  to  create  a  chart  of  flow  rate  versus 
pressure  differential  using  the  average  pressure  upstream  of  the  orifice  plate.  The 
average  pressure  differential,  combined  with  the  value  of  the  upstream  pressure  is  then 
used  to  estimate  the  average  flow  rate  from  the  chart.  Explanation  of  the  code  and  the 
code  itself  can  be  found  in  Appendix  C.  The  average  mass  flow  rate  of  the  fuel  is 
estimated  as  the  fuel  mass  lost  during  the  firing  divided  by  the  test  duration.  A  sample  of 
the  combustion  pressure  data  is  shown  in  Figure  17.  The  specific  impulse  is  determined 
by  dividing  the  total  impulse  by  the  weight  of  the  propellant  (both  fuel  and  oxidizer) 
consumed.  The  total  impulse  is  calculated  by  numerically  integrating  the  thrust  over  time 
for  the  duration  of  the  test  using  a  simple  Riemann  squares  approximation,  and  then 
divided  by  the  weight  of  the  propellant  used  as  shown  in  Eq.  13: 


2h^,+7;)*ft+a,-o 

(=o  2 


g  *  m 


prop 


(13) 


where  Isp  is  the  specific  impulse,  7)  is  the  current  thrust  value,  Ti+Al  is  the  thrust  at  the 
next  time  step,  t,  is  the  current  time,  ti+At  is  the  time  at  the  next  time  step,  g  is  the 

acceleration  due  to  gravity,  and  mprop  is  the  mass  of  the  propellant  consumed  during  the 
burn. 

The  major  uncertainty  in  the  specific  impulse  calculation  is  due  to  the  uncertainty 
in  the  oxidizer  mass  flow  rate  during  the  test.  This,  in  turn,  introduces  uncertainty  in  the 
mass  of  propellant  consumed.  The  burn  time  of  the  test  is  known  to  a  high  degree  of 
certainty,  as  noted  in  Table  8.  The  thrust  data  is  modified  to  remove  any  pressure  thrust 
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Combustion  Chamber  Pressure  vs  Time 
Test  173 


Figure  17:  Sample  Combustion  Chamber  Pressure  vs.  Time  Data 


contribution  to  remove  any  variations  in  the  overall  thrust  due  to  the  exit  pressure  not 
expanding  exactly  to  atmospheric  pressure.  An  example  of  the  thrust  data  collected  is 
shown  in  Figure  18.  This  is  accomplished  by  estimating  the  exit  pressure  at  every  point 
in  time  of  a  test  using  the  combustion  chamber  pressure  and  an  estimate  for  the  ratio  of 
the  exit  pressure  and  the  combustion  chamber  pressure  from  compressible  flow  dynamics 
as  shown  in  Eq.  14: 

pM)=—P'  W  (14) 

Pc 

where  pe(t)  is  the  exit  plane  pressure  as  a  function  of  time  and  pc(t)  is  the  combustion 
chamber  pressure  as  a  function  of  time.  Another  reason  for  this  thrust  correction  is  for 
comparison  to  the  analytical  data  produced.  The  ratio  of  exit  pressure  to  combustion 
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Thrust  vs  Time 

Beeswax  and  Oxygen,  Test  173 


Figure  18:  Sample  Thrust  vs.  Time  Data 


pressure  is  determined  from  the  nozzle  area  ratio  and  the  value  of  the  ratio  of  specific 
heats  (y).  The  nozzle  area  ratio  was  computed  using  Equation  15: 


A 
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(15) 


where  Ae  is  the  nozzle  exit  area,  A,  is  the  nozzle  throat  area,  de  is  the  nozzle  exit  diameter, 
and  dt  is  the  nozzle  throat  diameter.  The  nozzle  area  ratio  relation  is  then  solved  for  the 
exit  Mach  number  for  the  nozzle  as  shown  in  Equation  16: 
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where  Me  is  the  nozzle  exit  Mach  number.  Using  the  nozzle  exit  Mach  number  and  the 
ratio  of  specific  heats,  the  ratio  of  combustion  chamber  pressure  to  nozzle  exit  pressure 
can  be  calculated  as  shown  in  Equation  17: 


(17) 


The  value  of  the  ratio  of  specific  heats  was  determined  from  theoretical  analysis  of 
various  fuel  and  oxidizer  combinations  at  the  range  of  equivalence  ratios  used  in  the 
experimental  testing.  Table  9  shows  the  values  of  y  and  the  ratio  of  exit  pressure  to 
combustion  pressure  for  various  fucl/oxidizcr  combinations  for  a  combustion  pressure  of 
150  psig.  The  regression  rate  is  determined  over  a  range  of  equivalence  ratios.  The 
equivalence  ratio  (a)  is  a  measure  of  whether  the  fuc  l/oxidizcr  combination  is  rich,  lean, 
or  stoichiometric.  For  hybrid  rocket  motors,  it  is  defined  as  the  ratio  of  the  actual 
oxidizer  to  fuel  ratio  divided  by  the  stoichiometric  ratio  of  oxidizer  to  fuel  (Eq.  18). 


Given  that  the  combustion  occurs  in  a  diffusion  flame  and  the  fuel  and  oxidizer  are 


consumed  at  the  flame  front  in  the  boundary  layer  where  the  oxidizer  and  fuel  are  at  the 
stoichiometric  ratio,  the  equivalence  ratio  still  has  meaning.  Between  the  flame  front  and 
the  fuel  surface  is  a  fuel  rich  region  (31).  In  the  space  between  the  boundary  layer 
surface  and  the  flame  front  is  an  oxidizer  rich  region.  The  location  of  the  flame  front  in 
the  boundary  layer  is  determined  by  the  location  between  these  two  regions  where  the 
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Table  9:  Pressure  Ratios  for  Fuel/Oxidizer  Combinations 


Fuel/Oxidizer 

Y 

P  exit/P  combustion 

Lard/O? 

1.315 

0.082 

Beeswax/02 

1.306 

0.083 

Lard/N20 

1.306 

0.083 

Beeswax/N20 

1.304 

0.083 

equivalence  ratio  equals  one.  In  a  hybrid  rocket,  it  is  important  to  know  the  ratio  of  the 
oxidizer  input  to  the  motor  and  the  resulting  fuel  that  is  vaporized  from  the  grain  surface. 
This  allows  for  characterization  of  the  performance  of  the  hybrid  rocket  motor  as  a 
function  of  a  parameter  that  the  user  of  the  rocket  motor  has  physical  control  over.  The 
main  goal  of  the  regression  rate  analysis  is  to  detennine  the  regression  rate  as  a  function  a 
parameter  or  parameters  that  are  easily  controlled  or  measured  and  experimentally 
derived  constants.  Determining  the  regression  is  very  important  to  determining  the  thrust 
produced  by  the  motor.  This  is  due  to  the  mass  flow  rate  of  the  fuel  being  directly  related 
to  the  fuel  regression  rate.  Typical  parameters  used  are  total  mass  flux  rate,  axial 
distance  along  the  port,  length  of  the  port,  combustion  chamber  pressure,  oxidizer  mass 
flux  rate,  characteristic  exhaust  velocity,  and  oxidizer  mass  flow  rate  (32).  At  least  eight 
different  possible  relations  which  have  been  tried,  and  previous  studies  have  found  a  very 
low  degree  of  dependence  of  regression  rate  on  combustion  chamber  pressure.  However, 
it  has  been  shown  that  at  the  lower  end  of  combustion  pressures,  there  can  be  a  reduction 
in  the  regression  rate  due  to  the  radiation  heat  transfer  dominating  (1). 
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IV. V  Costs 


The  total  budget  for  this  research  was  approximately  $4500.  This  includes  the 
storage  building,  materials  to  construct  the  safety  wall,  instrumentation,  power  generator, 
beeswax,  lumber  to  construct  a  ramp  out  of  the  storage  building,  and  miscellaneous  parts 
(nuts,  bolts,  etc.).  This  budget  does  not  include  the  cost  of  building  the  test  stand  which 
was  constructed  during  the  1999-2000  school  year.  It  also  does  not  include  the  cost  of 
items  made  by  the  department  machine  shop  or  the  cost  of  filling  the  oxidizer  tanks. 
Table  10  below  shows  an  approximate  breakdown  of  the  costs. 


Table  10:  Research  Program  Expenses 


Item 

Cost 

Storage  Building 

$1500 

Power  Generator 

$500 

Safety  Wall  Materials 

$1000 

Beeswax 

$500 

Miscellaneous  Parts  (instrumentation, 
hardware,  electrical  equipment,  etc.) 

$1000 
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V.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


V.I  Repeatability  Study 

A  repeatability  study  was  performed  to  determine  how  consistently  the  test  stand 
and  instrumentation  performed.  It  is  important  to  determine  the  level  of  repeatability  in 
an  experimental  system  to  show  the  reliability  of  the  experimental  data.  A  high  level  of 
repeatability  improves  confidence  in  the  test  results.  If  poor  repeatability  is  observed, 
test  stand  problems  or  instrumentation  difficulties  could  be  indicated. 

For  this  study,  a  total  of  ten  repeatability  tests  were  performed  in  two  sets  of  five 
tests.  All  of  the  tests  were  run  with  beeswax  as  the  fuel  and  gaseous  oxygen  as  the 
oxidizer.  The  first  five  tests  used  an  oxidizer  supply  pressure  set  by  the  regulator  at  280 
psig.  The  oxidizer  supply  pressure  was  chosen  for  the  independent  variable  in  the 
repeatability  study  because  it  is  the  main  variable  in  the  experimental  setup  that  can  be 
controlled.  The  second  set  was  run  at  an  indicated  supply  pressure  of  500  psig.  These 
two  pressures  were  chosen  to  determine  repeatability  across  the  range  of  test  conditions 
used  in  this  study.  All  other  conditions  were  identical  in  all  of  the  tests.  The  time  for  the 
oxidizer  solenoid  was  set  at  5.0  s.  The  time  for  the  propane  solenoid/spark  plug  timer 
was  set  at  0.5  s.  The  thrust  and  combustion  chamber  pressure  were  recorded  and  then 
compared  for  each  set  of  tests.  Not  all  of  the  tests  were  run  on  the  same  day.  Table  1 1 
shows  the  ambient  temperature  on  the  day  for  each  of  the  ten  tests  run. 

A  summary  of  the  repeatability  tests  is  shown  in  Table  12.  For  the  tests  run  at  a 
supply  pressure  of  280  psig,  the  combustion  chamber  pressure  results  are  shown  in  Figure 
19.  Several  observations  can  be  made  from  Figure  20.  While  tests  123,  126,  132,  and 
136  are  fairly  close  in  chamber  pressure,  test  129  is  significantly  different. 
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Also,  test  126  shows  a  significant  spike  in  pressure  during  the  ignition  phase.  Overall, 
the  average  pressure  for  the  five  tests  was  approximately  102  psig  with  a  standard 
deviation  of  1 1  psig. 

The  thrust  data  for  the  five  tests  run  at  a  supply  pressure  of  280  psig  is  shown  in 
Figure  20.  The  first  observable  characteristic  is  the  significant  ignition  transients.  The 
first  0.5  s  of  each  test  show  a  large  spike  in  the  thrust  which  didn’t  begin  to  damp  out 
until  1.5  s.  Tests  129  and  132  were  close  in  value  and  shape.  These  two  tests  are  the  two 
steadiest  of  the  five  tests.  Test  126  has  an  oscillation  with  a  period  of  approximately  1.0 
s.  Test  123  has  an  oscillation  with  a  period  of  roughly  1.0  s.  and  an  amplitude  of 
approximately  35  N.  This  test  is  significantly  different  from  the  other  four  tests.  Test 
163  shows  a  general  decreasing  thrust  during  the  test,  but  has  a  thrust  value  roughly  10  N 
lower  than  test  129.  The  average  thrust  of  the  fives  tests  was  approximately  54  N  with  a 
standard  deviation  of  1 1  N. 


Combustion  Chamber  Pressure  vs  Time 
Repeatibility  Study,  Supply  Pressure  =  280  psi 


P 


Test  123 
Test  126 
Test  129 
Test  132 
Test  163 


Figure  19:  Comparison  of  Combustion  Pressures,  Supply  Pressure  =  280  psig 
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Thrust  vs  Time 

Repeatibility  Study,  Supply  Pressure  =  280  psig 


-Test  123 
-Test  126 
-Test  129 
-Test  132 
-Test  163 


Figure  20:  Comparison  of  Thrust,  Supply  Pressure  =  280  psig 


The  combustion  chamber  pressure  results  for  the  five  tests  run  with  a  supply 
pressure  of  500  psig  are  illustrated  in  Figure  2 1 .  All  five  tests  were  very  close  in  value  of 
the  pressure  and  trend.  The  average  combustion  chamber  pressure  of  the  lives  tests  was 
approximately  215  psig  with  a  standard  deviation  of  5  psig.  Overall,  the  combustion 
chamber  pressure  was  very  consistent  during  these  live  tests. 

Figure  22  shows  the  thrust  curve  comparison  of  the  live  repeatability  tests  run 
with  an  oxidizer  supply  pressure  of  500  psi.  Compared  to  the  combustion  pressure,  the 
thrust  curves  show  much  more  dispersion.  There  is  a  large  gap  between  tests  168,  170, 
trends,  but  the  values  of  the  thrust  vary  widely.  Test  170  has  an  increasing  thrust  versus 
and  171  and  tests  176  and  177.  Tests  168,  171,  176,  and  177  show  similar  constant  thrust 
time  trend,  which  is  significantly  different  from  the  other  four  tests.  The  average  thrust 
of  the  live  tests  was  approximately  140  N  and  the  standard  deviation  was  27  N.  Another 
characteristic  to  note  is  the  ignition  transients  and  the  time  it  took  for  them  to  damp  out. 


54 


Combustion  Chamber  Pressure 
(psig) 


Combustion  Chamber  Pressure  vs  Time 
Repeatibility  Study,  Supply  Pressure  =  500  psig 


- Test  168 

- Test  170 

- Test  171 

- Test  176 

- Test  177 


Time  (s) 


Figure  2 1 :  Comparison  of  Combustion  Pressures,  Supply  Pressure  =  500  psig 


Thrust  vs  Time 

Repeatibility  Study,  Supply  Pressure  =  500  psig 


- Test  168 

- Test  170 

- Test  171 

- Test  176 

- Test  177 


Figure  22:  Comparison  of  Thrust,  Supply  Pressure  =  500  psig 
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Tests  176  and  177  damped  out  the  quickest.  Tests  170  and  171  have  very  similar  ignition 
transients  despite  the  fact  that  test  170  has  an  increasing  thrust  and  test  171  has  a  constant 
thrust.  Test  168  had  the  largest  ignition  transients  of  the  five  tests.  However,  it  should  be 
noted  that  the  ignition  transients  are  much  less  severe  at  the  higher  supply  pressure  due  to 
much  more  stable  combustion. 

Due  to  the  oscillations  in  the  thrust  data,  an  attempt  was  made  to  determine  the 
natural  frequency  of  the  test  stand.  An  impulse  was  applied  to  the  test  stand  and  the 
response  was  measured.  There  were  three  such  tests.  The  time  response  data  from  the 
strain  gage  was  analyzed  using  a  fast  Fourier  transform  (FFT).  The  average  natural 
damped  frequency  of  the  three  tests  was  approximately  6.6  Hz. 

There  are  several  conclusions  from  these  results.  The  combustion  pressure 
becomes  more  repeatable  as  the  supply  pressure  increases.  The  problems  in  the  pressure 
data  with  the  tests  run  at  low  supply  pressures  could  be  due  to  combustion  instabilities 
caused  by  the  low  combustion  pressures.  The  spikes  in  tests  126  and  171  are  due  to 
inconsistencies  in  the  ignition  process.  This  shows  that  the  ignition  process  using 
propane  and  a  spark  plug  is  not  always  uniform.  Another  possible  cause  for  variances  in 
the  value  of  the  combustion  pressures  could  be  imprecise  setting  of  the  oxidizer  supply 
pressure.  It  is  quite  difficult  to  set  the  same  pressure  on  the  regulator  over  multiple  tests 
due  to  the  analog  pressure  indicator.  An  illustration  of  this  problem  is  seen  in  Figure  23. 
A  more  accurate  (i.e.  digital)  control  on  the  regulator  could  reduce  some  of  the 
differences  in  the  values  of  the  combustion  pressure.  As  with  the  combustion  pressure, 
the  thrust  becomes  more  stable  at  higher  pressures  due  to  more  stable  combustion  at  high 
pressure.  However,  there  is  still  a  wide  variance  in  the  values  of  the  thrust  even  at  the 
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Average  Combustion  Chamber  Pressure  vs  Oxidizer  Supply 

Pressure 

Beeswax  and  Oxygen 
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Figure  23:  Average  Combustion  Chamber  Pressure  vs.  Oxidizer  Supply  Pressure 
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higher  combustion  pressures.  It  is  apparent  that  due  to  the  inability  to  set  the  supply 
pressure  on  the  regulator,  this  is  not  the  best  way  to  examine  repeatability.  Therefore,  the 
supply  pressure  setting  can  only  be  used  to  choose  a  general  range  of  test  conditions  and 
the  combustion  chamber  pressures  have  to  be  examined  after  the  tests  to  determine  the 
actual  test  conditions.  A  better  way  to  examine  the  supply  pressure  is  to  plot  the  pressure 
upstream  of  the  orifice  plate  versus  time  (see  Figures  24  and  25).  The  average  supply 
pressure  in  Figure  24  is  approximately  1 10  psig  with  a  standard  deviation  of  9  psig.  The 
average  supply  pressure  in  Figure  25  is  approximately  230  psig  with  a  standard  deviation 
of  5  psig.  Another  way  to  examine  reliability  is  to  examine  the  entire  set  of  data  to  look 
for  general  trends.  One  option  is  to  examine  the  thrust  as  a  function  of  the  combustion 
chamber  pressure,  as  shown  in  Figure  26.  Aside  from  a  few  tests,  there  appears  to  be  a 
very  linear  trend  in  the  data.  This  would  indicate  that  the  test  stand  produces  more 
repeatable  results  than  shown  earlier.  A  strong  source  of  variance  in  the  results  is  the 
high  uncertainty  in  the  oxidizer  mass  flow  rate  calculation.  Other  sources  of  differences 
could  be  changes  in  ambient  temperatures,  uncertainty  in  the  fuel  mass  flow  rate,  and 
inconsistencies  (i.e  air  pockets,  possible  chemical  variations  between  batches  of  beeswax, 
etc.)  in  the  fuel  grains. 

V.II  Beeswax  and  Nitrous  Oxide 

As  stated  in  the  introduction  chapter,  one  of  the  oxidizers  chosen  for  examination 
was  nitrous  oxide.  It  is  a  commonly  used  oxidizer  and  is  very  safe  to  handle.  Several 
tests  were  run  in  an  attempt  to  acquire  useful  data.  However,  there  were  significant 
problems  with  the  testing.  None  of  the  tests  were  able  to  achieve  stable  combustion.  The 
worst  tests  failed  to  achieve  any  ignition.  The  best  tests  achieved  some  amount  of 
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Upstream  Pressure  vs  Time 
Repeatibility  Study,  Supply  Pressure  =  280  psig 
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Figure  24:  Upstream  Pressure  vs  Time,  Supply  Pressure  =  280  psig 


Upstream  Pressure  vs  Time 
Repeatibility  Study,  Supply  Pressure  =  500  psi 
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Figure  25:  Upstream  Pressure  vs  Time,  Supply  Pressure  =  500  psig 
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Average  Thrust  vs  Average  Combustion  Chamber  Pressure 
Beeswax  and  Oxygen 
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Figure  26:  Average  Thrust  vs  Average  Combustion  Chamber  Pressure 


ignition  and  had  significant  pulsing  during  the  burn.  The  exhaust  plumes  consisted 
mainly  of  smoke  with  little  visible  flame.  Figure  27  shows  some  representative  results  of 
the  testing.  Test  84  was  run  with  an  oxidizer  supply  pressure  of  240  psig  and  the  propane 
timer  was  set  at  1.0  s.  It  achieved  an  average  combustion  chamber  pressure  of 
approximately  55  psig,  but  did  not  achieve  stable  combustion.  At  the  time,  that  was  near 
the  limit  of  the  nitrous  oxide  regulator,  so  a  new  regulator  was  acquired  that  could  deliver 
a  supply  pressure  of  up  to  500  psig.  Test  92  was  run  at  a  supply  pressure  of  420  psig,  and 
the  propane  timer  was  increased  to  1.2  s.  The  propane  time  was  increased  from  1.0  s  to 
1 .2  s  to  provide  extra  initial  energy  to  the  fuel  grain  to  improve  the  ignition  process.  The 
average  pressure  was  approximately  35  psig.  Test  95  was  run  at  the  same  propane  time 
as  test  92,  but  the  supply  pressure  was  increased  to  500  psig.  The  average  combustion 
pressure  increased  to  only  40  psig.  Test  109  was  run  at  a  supply  pressure  of  500  psig,  but 
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Combustion  Chamber  Pressure  vs  Time 
Beeswax  and  Nitrous  Oxide 
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Figure  27:  Comparison  of  Combustion  Chamber  Pressures  for  Beeswax  and  Nitrous 


the  propane  time  was  increased  to  2.0  s.  This  produced  an  average  pressure  of 
approximately  60  psig.  The  length  of  this  test  was  increased  from  6.0  s  to  8.0  s  in  order 
to  give  the  ignition  process  longer  to  take  hold  and  initiate  the  combustion  process. 
However,  the  pressures  achieved  in  these  four  tests  were  too  low  for  stable  combustion  to 
occur.  Since  the  regulator  was  operating  at  its  maximum  discharge  pressure,  a  new 
solution  to  increase  the  combustion  chamber  pressure  had  to  be  formulated. 

An  investigation  of  the  characteristic  exhaust  velocity  was  initiated.  The 
characteristic  exhaust  velocity  is  a  measure  of  how  well  the  combustion  process  occurs  in 
a  rocket  motor  combustion  chamber.  It  is  experimentally  determined  and  shown  in  Eq. 

19: 

P  A 

c*  =  — - — —  (19) 

rii 
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where  c*  is  the  characteristic  exhaust  velocity,  Pc  is  the  combustion  pressure,  At  is  the 
nozzle  throat  area,  and  m  is  the  mass  flow  rate.  A  comparison  plot  of  the  theoretical 
characteristic  exhaust  velocity  for  beeswax/02  and  beeswax/N20  was  examined  and  is 
shown  in  Figure  28.  This  was  accomplished  using  the  Air  Force  Specific  Impulse  code, 
which  was  explained  in  the  Chapter  II.  The  combustion  chamber  pressure  examined  was 
150  psig,  and  the  exhaust  pressure  was  set  at  14.7  psig.  As  shown  in  Figure  28,  using 
nitrous  oxide  as  the  oxidizer  reduces  the  characteristic  exhaust  velocity  due  to  the 
presence  of  nitrogen  compared  to  using  pure  oxygen  as  the  oxidizer.  When  the  nozzle 
was  designed,  it  was  designed  assuming  that  oxygen  would  be  the  oxidizer.  If  the 
characteristic  exhaust  velocity  is  decreased  and  the  mass  flow  rate  is  assumed  constant, 
then  the  nozzle  throat  area  will  have  to  be  decreased  in  order  to  achieve  the  same 
combustion  chamber  pressure  as  shown  in  Eq.  20. 


Based  on  this,  a  new  nozzle  was  designed  for  use  with  nitrous  oxide.  Due  to  the 
restrictions  of  the  current  experimental  setup,  the  nozzle  had  to  be  the  same  length  as  the 
original  nozzle.  The  new  throat  diameter  was  0. 14  in.  and  the  new  exit  diameter  was  0.25 
in.  Test  112  was  run  with  a  supply  pressure  of  340  psig  and  the  propane  timer  set  at  2.0 
s.  The  result  was  a  maximum  combustion  chamber  pressure  of  150  psig.  This  was  a 
significant  increase  over  the  previous  four  tests.  However,  this  still  did  not  result  in 
stable  combustion,  and  the  motor  pulsed  during  the  test.  The  supply  pressure  was 
increased  to  380  psig  for  test  1 13,  and  the  maximum  combustion  chamber  pressure  was 


190  psig.  Again,  the  motor  pulsed  during  the  test  and  never  reached  steady  state 
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Comparison  of  Charactersitic  Exhaust  Velocity  vs 
Equivalence  Ratio,  Pc  =  150  psig 


Figure  28:  Comparison  of  Characteristic  Exhaust  Velocity  vs.  Equivalence  Ratio 


combustion.  Test  117  was  run  with  a  supply  pressure  of  420  psig,  and  the  maximum 
chamber  pressure  was  approximately  235  psig.  This  test  showed  significant  pulsing,  and 
the  chamber  pressure  was  much  less  smooth  than  tests  112  and  113.  Given  that  the  new 
regulator  was  again  nearing  its  maximum  supply  pressure,  the  next  logical  step  would 
have  been  to  design  a  new  nozzle  with  an  even  smaller  throat  diameter.  Unfortunately, 
the  combustion  chamber  pressure  was  beginning  to  approach  its  maximum  allowed  value. 
It  appears  that  a  combustion  pressure  of  higher  than  235  psig  is  required  for  stable 
combustion  to  be  realized.  The  current  test  stand  needs  to  be  modified  or  redesigned  to 
accommodate  higher  combustion  chamber  pressures  in  order  to  continue  this  course  of 
research.  Due  to  these  difficulties,  no  data  on  regression  rate  or  specific  impulse  was 
collected  for  beeswax  and  nitrous  oxide.  The  thrust  curves  for  the  tests  in  Figure  26  are 
shown  in  Figure  29.  From  the  appearance  of  the  exhaust  plume,  most,  it  not  all,  of  the 
thrust  produced  was  due  to  the  oxidizer  blowing  through  the  nozzle  and  the  burning  of 
the  propane  during  the  ignition  process.  It  is  unclear  why  tests  1 12,  1 13,  and  1 17  have  a 
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lower  thrust  level  than  tests  84,  92,  95,  and  109  even  though  the  chamber  pressure  was 
higher. 

V.III  Beeswax  and  Oxygen:  Regression  Rate 

For  the  regression  rate  analysis,  the  two  goals  were  to  compare  the  regression  rate 
of  beeswax  and  oxygen  with  other  hybrid  propellant  combinations  and  to  determine  an 
analytical  expression  for  the  regression  rate.  Two  of  the  propellant  combinations 
compared  are  HTPB/LOX  and  paraffin/oxygen;  Eqs.  21  and  22  show  their  expressions 
for  regression  rate  as  a  function  of  the  oxidizer  mass  flux  rate,  respectively  (12). 


r  =  0.146Got0681 

(21) 

r  =  0.488GOV°'62 

(22) 

Figure  30  shows  the  results  of  the  regression  rate  comparison.  Beeswax  and  oxygen  have 
a  significantly  higher  regression  rate  than  HTPB/LOX  over  the  range  of  oxidizer  mass 
flux  rates  examined.  Based  on  the  experimental  data,  the  regression  rate  for 
beeswax/oxygen  is  at  least  three  times  higher  than  HTPB/LOX.  This  was  similar  to 
results  found  in  previous  research  into  paraffin  and  oxygen  as  discussed  in  chapter  one. 
However,  it  was  surprising  to  discover  that  beeswax  and  oxygen  demonstrated  a  higher 
regression  rate  than  paraffin  and  oxygen.  The  beeswax  data  shows  a  wide  range  of 
spread,  so  it  is  difficult  to  determine  exactly  how  much  higher  the  regression  rate  is  over 
that  of  paraffin.  Ligure  30  shows  that  there  is  significant  uncertainty  in  the  mass  flux  rate 
due  to  the  uncertainty  in  the  oxidizer  mass  flow  rate.  The  uncertainty  in  the  regression 
rate  calculation  due  to  uncertainties  in  the  initial  port  radius,  burn  time,  initial  fuel  mass, 
and  burned  fuel  mass  was  investigated,  and  the  details  are  presented  in  Appendix  D. 
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Thrust  vs  Time 
Beeswax  and  Nitrous  Oxide 
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Figure  29:  Thrust  vs  Time,  Beeswax  and  Nitrous  Oxide 
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The  average  uncertainty  in  the  regression  rate  over  all  of  the  tests  was  0.33  mm/s.  It 
should  be  noted,  however,  that  despite  the  uncertainty  in  the  oxidizer  mass  flux  rate  and 
regression  rate,  the  regression  rate  is  still  significantly  higher  than  HTPB.  However,  the 
fact  that  the  regression  rate  is  higher  than  the  traditional  hybrid  propellants  such  as 
HTPB/LOX  shows  that  beeswax  should  be  able  to  provide  a  higher  thrust  hybrid  rocket 
motor  using  a  simple,  single  port  grain  geometry. 

The  next  step  was  to  determine  an  expression  for  the  regression  rate.  In  the 
literature,  there  are  many  different  expressions  for  regression  rate  as  a  function  of  many 
different  variables  (32).  Several  different  expressions  were  examined  to  evaluate  which 
would  best  fit  the  experimental  data.  One  thing  to  note  is  all  of  the  beeswax  fuel  grains 
tested  were  the  same  length,  so  the  length  of  the  combustion  port  was  not  considered  here 


Average  Regression  Rate  vs  Oxidizer  Mass  Flux  Rate 
Beeswax  and  Oxygen 
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Figure  30:  Comparison  of  Regression  Rates 


66 


as  a  dependent  variable  in  the  regression  rate  expressions.  This  is  due  to  the  nature  of  the 
current  test  stand  being  unable  to  accommodate  different  length  fuel  grains.  To  include 
the  effects  of  the  length  of  the  port,  different  length  fuel  grains  need  to  be  tested  and  a 
correlation  would  have  to  be  found.  For  this  study,  the  regression  rate  was  assumed 
constant  along  the  grain  length.  This  does  not  include  the  space  at  the  end  of  the  fuel 
grain  that  acts  as  the  aft  mixing  chamber.  All  of  the  expressions  examined  are  the 
typically  used  general  power  series,  shown  in  Eq.  23: 

r  =  aXn  (23) 

where  X  is  the  independent  variable  and  a  and  n  are  experimentally  derived  constants 
(1,2).  Four  independent  variables  were  analyzed  in  this  study:  oxidizer  mass  flux  rate, 
total  propellant  mass  flux  rate,  combustion  chamber  pressure,  and  oxidizer  mass  flow 
rate. 

As  shown  in  Eqs.  2 1  and  22,  it  is  common  to  examine  the  regression  rate  as  a 
function  of  the  oxidizer  mass  flux  rate.  A  power  series  expression  from  the  experimental 
data  in  Figure  30  is  shown  in  Eq.  24: 

r  =  1.573  1Got0'1026  (24) 

where  r  is  the  regression  rate  and  Gox  is  the  oxidizer  mass  flux  rate.  The  problem  with 
this  expression  is  the  poor  correlation  due  to  the  large  amount  of  scatter  in  the 
experimental  data.  The  R-squared  value  for  this  correlation  is  0.0054.  In  an  attempt  to 
improve  correlation  between  the  regression  rate  expression  and  the  experimental  data,  the 
total  propellant  mass  flux  rate  was  examined.  The  total  mass  flux  rate  includes  the  mass 
flux  of  the  fuel,  which  is  found  using  the  average  mass  flow  rate  of  the  fuel  over  the 
duration  of  the  test,  and  the  mass  flux  of  the  oxidizer.  The  resulting  expression  for  the 
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regression  rate  is  shown  in  Eq.  25 


r  =  0.0535Gtoi,14855  (25) 

where  Gtot  is  the  total  propellant  mass  flux  rate.  This  expression  has  a  much  better 
correlation  (R~  =  0.3642)  than  with  the  oxidizer  mass  flux  rate.  When  the  regression  rate 
is  plotted  against  the  total  propellant  mass  flux  rate,  as  in  Figure  31,  there  is  a  noticeable 
decrease  in  the  scatter  in  the  data.  While  the  total  propellant  mass  flux  rate  shows  a 
better  correlation  than  the  oxidizer  mass  flux  rate,  there  are  other  expressions  that  have 
better  agreement  with  the  experimental  data.  A  common  independent  variable  used  for 
regression  rate  expressions  in  solid  rocket  motors  is  the  combustion  chamber  pressure. 
The  expression  for  regression  rate  as  a  function  of  the  combustion  chamber  pressure  is 
shown  in  Eq.  26 

r  =  0.18LP,0- 478  (26) 

where  Pc  is  the  combustion  chamber  pressure.  The  reduction  in  the  scatter  of  the 
regression  rate  data  when  presented  against  the  combustion  chamber  pressure  can  be  seen 
in  Figure  32.  The  R-squared  value  for  this  correlation  is  0.4758.  The  expression  in  Eq. 

12  was  evaluated  only  using  the  data  shown  in  Figure  32  that  used  a  converging- 
diverging  nozzle.  The  data  resulting  from  using  a  converging  nozzle  is  shown  for 
historical  reference  only. 

The  last  parameter  examined  as  an  independent  variable  was  the  oxidizer  mass 
flow  rate.  Figure  33  shows  the  regression  rate  versus  the  oxidizer  mass  flow  rate.  Using 
the  oxidizer  mass  flow  rate  showed  the  least  scatter  of  all  of  the  variables  analyzed. 

When  an  expression  for  the  regression  rate  was  determined,  it  showed  the  best  correlation 
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Regression  Rate  vs  Total  Propellant  Mass  Flux  Rate 
Beeswax  and  Oxygen 


Total  Propellant  Mass  Flux  Rate  (g/cm2  s) 


Figure  3 1 :  Regression  Rate  vs.  Total  Propellant  Mass  Flux  Rate 


Regression  Rate  vs  Combustion  Chamber  Pressure 
Beeswax  and  Oxygen 


Figure  32:  Regression  Rate  vs.  Combustion  Chamber  Pressure 


69 


with  the  experimental  data  (R2  =  0.4876).  The  analytical  expression  is  shown  in  Eq.  27: 

r  =  1 1.815mox0'4862  (27) 

where  mox  is  the  oxidizer  mass  flow  rate. 

While  using  the  oxidizer  mass  flow  rate  has  the  best  correlation,  there  is  still  significant 
scatter  in  the  data.  The  maximum  spread  in  the  regression  rate  at  a  given  mass  flow  rate 
is  ±  0.4  mm/s. 

There  are  several  possible  reasons  for  the  significant  scatter  in  the  regression  rate 
data.  As  stated  in  Chapter  IV,  the  error  from  the  pressure  transducers  was  ±  0.4%.  While 
this  is  sufficient  for  the  measurement  of  the  pressure  upstream  and  downstream  of  the 
orifice  plate,  the  mass  flow  rate  is  computed  from  the  pressure  differential  across  the 
orifice  plate.  When  measuring  the  difference  in  pressure,  a  0.4%  error  in  the  transducers 
begins  to  have  a  much  more  significant  effect  on  the  data.  This  is  illustrated  in  Figure  34. 


Regression  Rate  vs  Oxidizer  Mass  Flow  Rate 
Beeswax  and  Oxygen 


Figure  33:  Regression  Rate  vs.  Oxidizer  Mass  Flow  Rate 
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The  pressure  differential  across  the  orifice  plate  is  shown  versus  time  for  Test 
169.  This  test  had  an  oxidizer  supply  pressure  of  400  psig.  Another  possible  causes  for 
scatter  in  the  data  include  inconsistencies  in  the  fuel  grains  due  to  air  pockets  formed 
during  the  fabrication  process.  One  way  to  measure  this  possibility  is  to  examine  the 
initial  fuel  grain  mass  measurements.  The  average  initial  fuel  grain  mass  was  3.391  kg 
with  a  standard  deviation  of  0.017  kg.  Another  cause  could  be  the  spallation  of  solid  or 
melted  particles  of  beeswax  leaving  the  motor  unburned.  This  would  tend  to  inflate  the 
measured  regression  rate,  as  no  attempt  to  measure  the  spallation  of  unburned  beeswax 
was  made.  Another  possibility  is  errors  in  the  fuel  grain  mass  measurements.  As  noted 
earlier,  the  accuracy  of  the  scale  used  was  ±  lgram.  All  of  the  tests  were  run  for  the  same 
duration  (5.0  s)  and  the  propane  was  run  for  the  same  duration  (0.5  s)  for  all  tests.  The 
likely  explanation  of  the  errors  in  the  data  is  a  combination  of  all  of  these  possibilities. 


Pressure  Difference  Across  Orifice  Plate  vs  Time 
Beeswax  and  Oxygen,  Test  169 
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Figure  34:  Pressure  Difference  vs.  Time,  Test  169 
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V.IV  Beeswax  and  Oxygen:  Specific  Impulse 

The  sea  level  specific  impulse  for  beeswax  and  oxygen  was  found  by  integrating 
the  thrust  over  the  duration  of  the  test  to  find  the  total  impulse  and  then  dividing  by  the 
weight  of  the  propellant  burned.  The  experimental  data  was  used  to  evaluate  the  specific 
impulse,  and  then  the  data  was  compared  to  a  theoretical  analysis  of  the  specific  impulse. 
The  theoretical  data  was  evaluated  at  a  combustion  chamber  pressure  of  150  psig,  and  the 
equivalence  ratio  was  varied.  The  theoretical  seal  level  specific  impulse  was  computed 
assuming  frozen  flow  through  the  nozzle.  A  plot  of  the  experimental  and  theoretical 
specific  impulses  as  a  function  of  equivalence  ratio  is  presented  in  Figure  35.  Several 
observations  can  be  made  from  Figure  35.  First,  a  significant  variance  in  the  data  is 
apparent.  The  primary  source  of  scatter  is  the  uncertainty  in  the  equivalence  ratio,  which 
is  primarily  due  to  uncertainty  in  the  oxidizer  mass  flow  rate.  A  calculation  of  the 
uncertainty  in  the  equivalence  ratio  was  performed  and  the  average  uncertainty  over  all 
the  tests  was  0.24.  This  is  illustrated  in  Figure  35  by  horizontal  error  bars.  This  again 
highlights  the  need  for  funding  to  purchase  a  differential  pressure  transducer  to  improve 
accuracy  of  the  oxidizer  mass  flow  rate.  Also,  a  large  number  of  the  tests  resulted  in  a 
calculated  specific  impulse  greater  than  the  theoretical  value.  There  are  several  possible 
reasons  for  these  results.  Clearly,  there  is  enough  scatter  in  the  data  to  question  the 
specific  impulse  data  that  is  above  the  theoretical  curve  in  Figure  35.  The  possible 
sources  of  uncertainty,  in  descending  order  of  magnitude,  are  the  oxidizer  mass  flow  rate, 
the  burned  fuel  mass,  the  burn  time,  and  the  thrust.  When  all  of  these  are 
accounted  for,  the  average  uncertainty  over  all  the  tests  for  the  calculated  sea  level 
accounted  for,  the  average  uncertainty  over  all  the  tests  for  the  calculated  sea  level 
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specific  impulse  was  5.6  s.  The  calculations  for  the  specific  impulse  uncertainty  are 
shown  in  Appendix  D.  Another  variable  is  the  combustion  pressure.  The  experimental 
data  was  collected  over  a  range  of  combustion  chamber  pressures,  while  the  theoretical 
analysis  was  conducted  for  a  chamber  pressure  of  150  psig.  However,  if  the  specific 
impulse  is  analyzed  as  a  function  of  the  combustion  chamber  pressure,  similar  results  are 
observed  as  shown  in  Figure  36.  There  are  several  reasons  one  would  expect  the 
experimental  specific  impulse  data  to  be  below  the  calculated  theoretical  values.  The 
theoretical  values  do  not  account  for  a  chemically  reacting  gas  throughout  the  nozzle.  In 
other  words,  the  computer  code  assumes  the  chemical  composition  of  the  gas  does  not 
change  as  it  flows  through  the  nozzle.  This  can  change  the  ratio  of  specific  heats  and  the 
molecular  mass  of  the  exhaust  gas  which  in  turn  changes  the  exit  velocity.  This  can 
represent  a  loss  of  specific  impulse  of  0.5%  (2).  Also  not  accounted  for  in  the  theoretical 
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Figure  35:  Comparison  of  Experimental  and  Theoretical  Specific  Impulse 
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analysis  are  viscous  losses  in  the  nozzle,  which  reduce  the  exit  velocity  and  in  turn  reduce 
the  specific  impulse.  These  losses  are  due  to  the  presence  of  a  boundary  layer  inside  the 
nozzle,  and  can  reduce  the  effective  exhaust  velocity  by  0.5%  to  1.5%  (1).  Multiphase 
flow  in  the  nozzle  can  also  reduce  the  exit  velocity  as  discussed  previously.  Having 
liquid  or  solid  particles  in  the  exhaust  can  cause  a  reduction  of  up  to  5%  in  the  specific 
impulse  of  the  motor.  Other  losses  can  be  caused  by  unsteady  or  oscillating  combustion. 
As  shown  in  the  repeatability  study  earlier,  there  was  unsteady  combustion  occurring  in 
some  of  the  tests,  primarily  at  lower  combustion  pressures  (2).  A  final  source  of 
reduction  in  the  specific  impulse  is  divergence  losses  due  to  the  shape  of  the  nozzle.  Due 
to  the  conical  shape  of  the  nozzle,  there  will  be  a  component  of  the  exit  velocity  vector  in 
the  radial  direction.  Any  component  of  the  exit 
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Figure  36:  Specific  Impulse  vs.  Combustion  Chamber  Pressure 
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velocity  vector  not  in  the  axial  direction  will  result  in  a  loss  of  specific  impulse.  This 
contribution  to  the  losses  is  a  function  of  the  cosine  of  the  nozzle  half  angle,  and  can  be 
computed  using  Eq  28  (1): 

2  =  t(l  +  cos0j  (28) 

where  X  is  the  thrust  efficiency  and  0cn  is  the  nozzle  half  angle.  The  half  angle  for  the 
nozzle  used  was  3.54  degrees,  which  results  in  a  loss  of  approximately  0.1%.  Adding  up 
these  losses  results  in  a  possible  7.1%  reduction  in  the  specific  impulse  from  the 
theoretical  values.  Figure  37  shows  these  losses  in  relation  to  the  theoretical  and 
experimentally  derived  results  for  the  specific  impulse.  Note  that  this  new  curve  for  the 
specific  impulse  does  not  account  for  possible  errors  in  the  propellant  mass  flow  rate  or 
thrust  measurements.  Ideally,  the  experimental  data  would  fall  very  close  to,  or  on,  the 
curve  that  accounts  for  the  losses.  There  is  also  the  possibility  of  unaccounted  losses  that 
would  be  responsible  for  the  discrepancies  between  the  experimental  data  and  the  loss- 
accounted  theoretical  results. 

V.V  Beeswax  and  Oxygen:  Thrust 

An  analysis  of  the  thrust  data  was  conducted  to  determine  a  theoretical  value  for 
the  thrust  for  a  hybrid  rocket  motor  using  beeswax  and  oxygen  as  propellants.  For 
simplicity,  it  was  assumed  the  combustion  chamber  pressure  was  150  psig  and  the 
equivalence  ratio  was  0.3.  This  combustion  pressure  was  chosen  because  the  part  of  the 
theoretical  analysis  used  to  compare  to  the  experimental  data  was  performed  at  that 
pressure.  This  equivalence  ratio  was  chosen  as  an  average  value  experienced  in  the 
experimental  results.  For  this  analysis,  it  was  also  assumed  to  have  an  ideally  expanded 
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nozzle,  which  eliminated  the  exit  pressure  contribution  to  the  thrust.  This  reduces  the 
thrust  expression  to  Eq.  29: 


T  =  mV e  (29) 

where  T  is  the  thrust,  m  is  the  mass  flow  rate  through  the  nozzle,  and  Ve  is  the  nozzle 
exit  velocity.  Using  an  ideally  expanded  nozzle  also  allows  the  use  of  the  relation  of  the 
specific  impulse  to  the  exit  velocity  as  shown  in  Eq.  30: 

q,=—  (30) 

go 

where  Isp  is  the  specific  impulse,  Ve  is  the  nozzle  exit  velocity,  and  go  is  acceleration  due 
to  gravity.  Substituting  Eq.  30  into  Eq.  29  yields  an  expression  for  the  thrust  as  shown  in 
Eq.  31. 
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Figure  37:  Specific  Impulse  Including  Losses 
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T  =  ™ISpg  0 


(31) 


Eq.  3 1  was  used  to  determine  the  theoretical  values  of  thrust.  From  Appendix  A,  the 
specific  impulse  for  beeswax/Cf  at  an  equivalence  ratio  of  0.3  is  184  s.  To  correlate  with 
the  experimental  data,  the  total  approximated  propellant  mass  flow  rate  from  the 
experimental  data  was  used  in  Eq.  3  E  A  comparison  of  the  measured  average  thrust 
values  with  their  theoretical  counterparts  is  shown  in  Table  13.  It  is  evident  from  Table 
13  that  the  experimental  values  of  the  thrust  are  significantly  higher  than  the  theoretical 
estimations.  This  is  likely  due  to  the  uncertainty  in  the  propellant  mass  flow  rates  used  in 
the  ideal  thrust  calculations.  It  would  appear  the  true  mass  flow  rates  are  higher  than 
have  been  approximated  since  the  actual  measured  thrust  levels  are  higher  than  the 
calculated  ideal  thrust  levels.  Higher  actual  values  of  equivalence  ratio  would  increase 
both  specific  impulse  and  propellant  mass  flow  rates  used  in  the  thrust  calculations  and 
more  closely  match  them  with  the  measured  thrust  values. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


VI.I  Conclusions 

The  main  objective  of  this  research  was  to  evaluate  regression  rates  for  non- 
conventional  bio-derived  hybrid  rocket  propellants  and  compare  them  to  traditional 
hybrid  rocket  motor  propellants.  Such  preliminary  work  is  necessary  to  produce  a  hybrid 
rocket  with  a  high  thrust  level  using  a  simple,  single  port  geometry.  For  the  tests  using 
beeswax  and  oxygen,  the  regression  rate  was  measured  and  determined  to  be  at  least  3 
times  higher  than  HTPB/LOX  and  possibly  higher  than  paraffin/oxygen.  For  example,  at 
an  oxidizer  mass  flux  rate  of  4.0  g/cm2s,  HTPB/LOX  has  a  regression  rate  of  0.4  min/s, 
paraffin/oxygen  has  a  regression  rate  of  1.2  min/s,  and  beeswax/oxygen  has  a  regression 
rate  of  1 .8  ±  0.33  mm/s.  The  wide  range  given  for  the  beeswax/oxygen  regression  rate  is 
due  to  scatter  in  the  experimental  data.  An  expression  for  the  regression  rate  was  derived 
as  a  function  of  the  oxidizer  mass  flow  rate.  Mass  flow  rate  was  chosen  as  the 
independent  variable  for  the  expression  because  it  had  the  best  correlation  with  the 
experimental  data.  Unfortunately  this  limits  the  application  of  the  equation  to  our  engine 
design  or  one  of  identical  design. 

Another  objective  was  to  determine  experimentally  the  thrust  and  specific  impulse 
and  compare  them  to  theoretical  values  to  validate  the  data.  This  part  of  the  research 
proved  difficult  due  to  mechanical  issues  resulting  in  poor  thrust  data.  This  led  to 
significant  scatter  in  the  thrust  data,  which  in  turn  produced  scatter  in  the  specific  impulse 
data.  Calculated  experimental  sea  level  specific  impulses  were  both  above  and  below 
their  theoretical  levels.  Some  tests  did  fall  near  the  theoretical  values  when  losses  were 
accounted  for,  but  many  tests  were  above  the  theoretical  value.  The  thrust  data  collection 
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needs  to  be  improved  in  order  to  reduce  uncertainty.  At  a  combustion  chamber  pressure 
of  150  psig,  the  experimentally-measured  thrust  turned  out  to  be  higher  than  the  predicted 
values.  The  test  stand  appears  to  be  reliable  for  combustion  pressure  measurements,  but 
improvements  are  needed  in  the  thrust  data  system.  The  largest  room  for  improvement  is 
in  the  measurement  of  the  oxidizer  mass  flow  rate.  Due  to  the  significant  uncertainties  in 
the  oxidizer  mass  flow  rate,  the  equivalence  ratios  calculated  in  this  study  have  a  very 
significant  uncertainty  and  should  be  viewed  accordingly. 

When  beeswax  and  nitrous  oxide  were  tested,  it  was  determined  that  a  higher 
combustion  chamber  pressure  than  the  test  stand  can  achieve  is  required  for  stable 
combustion.  Attempts  were  made  to  increase  the  combustion  chamber  pressure  by 
redesigning  the  nozzle,  but  the  pressure  limits  of  the  test  stand  were  reached  before  stable 
combustion  occurred.  No  data  for  regression  rate  or  specific  impulse  was  acquired,  but 
this  propellant  combination  is  still  worth  pursuing.  A  summary  of  all  the  tests  is 
presented  in  Appendix  F. 

For  the  analytical  model,  the  goals  were  to  produce  theoretical  data  at  a  chamber 
pressure  of  150  psig  for  comparison  with  the  experimental  data  and  at  a  chamber  pressure 
of  500.38  psig  for  comparison  with  published  data  for  traditional  hybrid  propellants.  A 
set  of  data  at  a  pressure  of  150  psig  was  created  and  used  to  validate  the  experimental 
data.  A  set  of  data  for  many  different  non-traditional  hybrid  propellants  was  produced 
for  a  pressure  of  500.38  psig  and  compared  to  traditional  hybrid  propellants.  Beeswax 
combined  with  oxygen  proved  to  be  the  best  propellant  combination  producing  a 
maximum  specific  impulse  of  327  s.  Overall,  oxygen  as  the  oxidizer  produced  better 
performance  than  nitrous  oxide.  For  the  solid  fuels  examined,  beeswax  produced  the  best 
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results  followed  by  lard.  After  lard  was  the  50%/50%  mixture  of  lard  and  paraffin  by 
mass.  When  comparing  beeswax  and  oxygen  with  the  traditional  hybrid  propellants,  it 
had  a  higher  specific  impulse  than  HTPB/LOX,  HTPB/hydrogen  peroxide,  and 
HTPB/nitrogen  tetroxide.  A  database  of  propellant  properties  was  created  for  use  by 
future  hybrid  rocket  motor  designers. 

Overall,  the  state  of  the  art  in  hybrid  rocket  motor  technology  was  advanced  by 
examining  previously  not  considered  non-conventional,  bio-derived  propellants. 

VI.II  Recommendations  for  Further  Research 

For  future  research,  several  improvements  to  the  test  stand  are  recommended. 

The  main  recommendation  is  to  improve  the  oxidizer  mass  flow  rate  measurement.  For 
the  pressure  data  used  to  determine  the  oxidizer  mass  flow  rate,  a  differential  pressure 
transducer  is  recommended  to  reduce  the  scatter  in  the  mass  flow  rate  measurements. 
Unfortunately,  the  funding  was  not  available  to  purchase  the  needed  differential  pressure 
transducer.  This  would  improve  the  correlation  for  the  regression  rate  expressions.  If  a 
differential  pressure  transducer  is  not  successful,  a  flow  rate  meter  could  be  acquired. 

The  effects  of  ignition  transients  and  oscillations  need  to  be  minimized  with  regard  to 
thrust  measurements.  A  preload  could  be  placed  on  the  system  during  testing.  This 
could  potentially  help  with  the  scatter  in  the  data  and  the  damping  characteristics.  If  this 
is  not  sufficient  to  correct  the  scatter  in  the  data,  then  the  current  strain  gage  setup  could 
be  replaced  with  a  load  cell.  Also,  the  test  stand  needs  to  be  able  to  accommodate  higher 
combustion  pressures.  The  nozzle  end  cap  could  be  modified  to  incorporate  an  o-ring  to 
assist  in  the  containing  the  pressure.  Also,  another  method  for  securing  the  fuel  grain  in 
the  test  stand  needs  to  be  developed.  The  practice  of  using  bolts  to  secure  the  fuel  grain 
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can  only  be  improved  by  using  stronger  bolts.  This  has  limits  as  the  current  setup  will 
not  handle  bolts  much  larger  than  currently  used.  Another  parameter  desirable  to 
measure  is  the  combustion  chamber  temperature  and  compare  it  to  the  theoretical  value. 
An  internet  search  was  unable  to  find  a  thermocouple  that  could  handle  the  maximum 
temperature  of  roughly  3400K.  However,  there  are  thermocouples  available  that  can 
measure  the  nozzle  exit  plane  temperature,  and  the  exit  temperature  can  be  used  to 

estimate  the  combustion  chamber  temperature  using  Equation _ .  Unfortunately,  these 

thermocouples  are  well  outside  the  funding  available  for  this  project.  A  lesser 
thermocouple  was  purchased  and  an  attempt  was  made  to  measure  the  exhaust 
temperature.  The  thermocouple  did  not  collect  any  data  due  to  the  strong  oxidizing 
environment,  but  it  physically  survived  the  temperature  and  dynamic  pressure.  A  list  of 
recommended  equipment  is  included  in  Appendix  E. 

For  future  research,  more  tests  need  to  be  performed  to  get  better  correlation  for 
the  regression  rate  data  and  thrust  data.  Improving  the  thrust  data  will  improve  the 
resulting  specific  impulse  data.  Also,  the  current  research  was  limited  to  a  range  of 
equivalence  ratios  between  0.2  and  0.4.  The  mass  flow  rate  of  the  oxidizer  system  needs 
to  be  increased  in  order  to  increase  the  equivalence  ratio.  It  is  desirable  to  collect  data 
over  a  wider  range  of  equivalence  ratios.  Once  these  tests  are  completed,  scale-up  hybrid 
motors  will  need  to  be  tested  in  order  to  determine  if  the  results  acquired  for  lab-scale 
motors  will  translate  into  flight-sized  hybrid  rocket  motors.  One  way  to  assist  in 
applying  the  lab  scale  tests  to  larger  motors  is  to  incorporate  the  length  of  the  combustion 
port  into  the  regression  rate  expression.  To  accomplish  this,  fuel  grains  of  differing 
lengths  need  to  be  tested.  Another  area  that  could  be  pursued  in  the  experimental  phase 
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is  the  use  of  additives  to  the  solid  fuels.  Aluminum  powder  could  be  added  to  evaluate  if 
it  matches  the  theoretical  performance.  Another  additive  that  could  be  used  is  carbon 
black.  This  is  often  added  to  solid  propellants  to  reduce  radiation  heat  transfer  to  the 
solid  fuel  grain. 
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APPENDIX  A:  DATA  FOR  CHAMBER  PRESSURE  OF  150  PSIG 


Temperature  vs  Equivalence  Ratio:  Beeswax  and  Oxygen 
(Pc  =  150  psi) 


Figure  38:  Temperature  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  150  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  150  psi) 


- Vac  Isp  (sec)  FF 

- Sea  Lvl  Isp  (sec)  FF 


Figure  39:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  150  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =150  psi) 


Figure  40:  Ratio  of  Specific  Fleats  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  150 
psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  150  psi) 


Figure  41:  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  150  psig 
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Temperature  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  42:  Temperature  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc  =  150  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =150  psi) 


- Vac  Isp  (sec)  FF 

- Sea  Lvl  Isp  (sec)  FF 


Figure  43:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc  =  150 
psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  44:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc 
150  psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  45:  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc  =  150 
psig 
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Temperature  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =  150  psi) 


Figure  46:  Temperature  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc  =  150  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =150  psi) 


Figure  47:  Specific  Impulse  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc  =  150  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =  150  psi) 


Figure  48:  Ratio  of  Specific  Fleats  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc  =  150  psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =  150  psi) 


Figure  49:  Molecular  Mass  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc  =  150  psig 
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Temperature  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  50:  Temperature  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  =  150  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  150  psi) 


Equivalence  Ratio 


- Vac  Isp  (sec)  FF 

[  Sea  Lvl  Isp  (sec)  FF 


Figure  51:  Specific  Impulse  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  =  150  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  52:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  =  150 
psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  53:  Molecular  Mass  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  =  150  psig 


97 


Temperature  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =150  psi) 


Figure  54:  Temperature  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc=  150  psig 


Specific  Impulse  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =150  psi) 


- Vac  Isp  (sec)  FF 

Sea  Lvl  Isp  (sec)  FF 


Figure  55:  Specific  Impulse  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc  =  150  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =  150  psi) 


Figure  56:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc  =  150 
psig 


Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =  150  psi) 


Figure  57:  Molecular  Mass  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc  =  150  psig 
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Temperature  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  58:  Temperature  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc  =  150  psig 

Specific  Impulse  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =150  psi) 


- Vac  Isp  (sec)  FF 

- Sea  Lvl  Isp  (sec)  FF 


Figure  59:  Specific  Impulse  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc  =  150 
psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =150  psi) 


Figure  60:  Ratio  of  Specific  Fleats  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc 
150  psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =  150  psi) 


Figure  61:  Molecular  Mass  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc  =  150 
psig 
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APPENDIX  B:  DATA  FOR  CHAMBER  PRESSURE  OF  500.38  PSIG 


Temperature  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  500.38  psi) 


Figure  62:  Temperature  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  500.38  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  500.38  psi) 


- Vac  Isp  (sec)  FF 

- Sea  Lvl  Isp  (sec)  FF 


Figure  63:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  500.38  psi) 


Figure  64:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  500.38 
psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  and  Oxygen  (Pc  =  500.38  psi) 


Figure  65:  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax/Oxygen,  Pc  =  500.38  psig 
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Temperature  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  66:  Temperature  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc  =  500.38  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  67:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc  =  500.38 
psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  68:  Ratio  of  Specific  Fleats  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc  = 
500.38  psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  69:  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax/Nitrous  Oxide,  Pc  =  500.38 
psig 
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Specific  Impulse  (sec) 


Temperature  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =  500.38  psi) 


Equivalence  Ratio 


Figure  70:  Temperature  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc  =  500.38  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =  500.38  psi) 


Equivalence  Ratio 


- Vac  Isp  (sec)  FF 

- Sea  Lvl  Isp  (sec)  FF 


Figure  71:  Specific  Impulse  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc=  500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =  500.38  psi) 


Equivalence  Ratio 


Figure  72:  Ratio  of  Specific  Fleats  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc  =  500.38  psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Lard  and  Oxygen  (Pc  =  500.38  psi) 


Figure  73:  Molecular  Mass  vs.  Equivalence  Ratio,  Lard/Oxygen,  Pc  =  500.38  psig 
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Temperature  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  74:  Temperature  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  =  500.38  psig 

Specific  Impulse  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


- Vac  Isp  (sec)  FF 

- Sea  Lvl  Isp  (sec)  FF 


Figure  75:  Specific  Impulse  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  =  500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  76:  Ratio  of  Specific  Fleats  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  = 
500.38  psig 


Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Lard  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  77:  Molecular  Mass  vs.  Equivalence  Ratio,  Lard/Nitrous  Oxide,  Pc  =  500.38  psig 
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Temperature  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =  500.38  psi) 


Figure  78:  Temperature  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc  =  500.38  psig 

Specific  Impulse  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =  500.38  psi) 


■Vaclsp(sec)  FF 
“Sea  Lvl  Isp  (sec)  FF  | 


Figure  79:  Specific  Impulse  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc  =  500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =  500.38  psi) 


Figure  80:  Ratio  of  Specific  Fleats  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc  =  500.38 
psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Oxygen  (Pc  =  500.38  psi) 


Figure  81:  Molecular  Mass  vs.  Equivalence  Ratio,  50/50  LP/Oxygen,  Pc  =  500.38  psig 
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Temperature  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  82:  Temperature  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc  =  500.38  psig 

Specific  Impulse  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


- Vac  Isp  (sec)  FF 

- Sea  Lvl  Isp  (sec)  FF 


Figure  83:  Specific  Impulse  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc  =  500.38 
psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  84:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc  = 
500.38  psig 

Combustion  Products  Molecular  Weight  vs  Equivalence  Ratio 
50%  Lard/50%  Paraffin  and  Nitrous  Oxide  (Pc  =  500.38  psi) 


Figure  85:  Molecular  Mass  vs.  Equivalence  Ratio,  50/50  LP/Nitrous  Oxide,  Pc  =  500.38 
psig 
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Temperature  vs  Equivalence  Ratio 
Beeswax  +  5%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  86:  Temperature  vs.  Equivalence  Ratio,  Beeswax  +  5%  Al/Oxygcn,  Pc  =  500.38 
psig 

Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  +  5%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  87:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax  +  5%  Al/Oxygcn,  Pc  = 
500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  +  5%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  88:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Beeswax  +  5%  Al/Oxygcn,  Pc 
=  500.38  psig 


Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  +  5%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  89:  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax  +  5%  Al/Oxygcn,  Pc  = 
500.38  psig 
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Temperature  vs  Equivalence  Ratio 
Beeswax  +  10%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  90:  Temperature  vs.  Equivalence  Ratio,  Beeswax  +  10%  Al/Oxygcn,  Pc=  500.38 
psig 


Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  +  10%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  91:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax  +  10%  Al/Oxygcn,  Pc  = 
500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  +  10%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  92:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Beeswax  +  10%  Al/Oxygcn,  Pc 
=  500.38  psig 


Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  +  10%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  93:  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax  +  10%  Al/Oxygcn,  Pc  = 
500.38  psig 
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Temperature  vs  Equivalence  Ratio 
Beeswax  +  15%  Al  and  Oxygen,  Pc  =  500.38  psi 


Equivalence  Ratio 

Figure  94:  Temperature  vs.  Equivalence  Ratio,  Beeswax  +  15%  Al/Oxygcn,  Pc  =  500.38 
psig 

Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  +  15%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  95:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax  +  15%  AEOxygen,  Pc  = 
500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  +  15%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  96:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Beeswax  +  15%  Al/Oxygcn,  Pc 
=  500.38  psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  +  15%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  97:  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax  +  15%  Al/Oxygcn,  Pc  = 
500.38  psig 
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Temperature  vs  Equivalence  Ratio 
Beeswax  +  20%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  98:  Temperature  vs.  Equivalence  Ratio,  Beeswax  +  20%  Al/Oxygcn,  Pc  =  500.38 
psig 

Specific  Impulse  vs  Equivalence  Ratio 
Beeswax  +  20%  Al  and  Oxygen,  Pc  =  500.38  psi 


Figure  99:  Specific  Impulse  vs.  Equivalence  Ratio,  Beeswax  +  20%  Al/Oxygcn,  Pc  = 
500.38  psig 
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Ratio  of  Specific  Heats  vs  Equivalence  Ratio 
Beeswax  +  20%  Al  and  Oxygen,  Pc  =  500.38  psi 


Equivalence  Ratio 

Figure  100:  Ratio  of  Specific  Heats  vs.  Equivalence  Ratio,  Beeswax  +  20%  Al/Oxygcn, 
Pc=  500.38  psig 

Combustion  Products  Molecular  Mass  vs  Equivalence  Ratio 
Beeswax  +  20%  Al  and  Oxygen,  Pc  =  500.38  psi 


Equivalence  Ratio 

Figure  101 :  Molecular  Mass  vs.  Equivalence  Ratio,  Beeswax  +  20%  AEOxygen,  Pc  = 
500.38  psig 
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APPENDIX  C:  MASS  FLOW  RATE  CALCULATIONS 


Equations  used  in  oxidizer  mass  flow  rate  chart  calculation: 

The  mass  flow  rate  is  calculated  as: 
e*C*  At*  yj(2*  p*  AP) 


m  = 


('-?) 


(32) 


Beta  ratio  is  the  ratio  of  orifice  plate  hole  diameter  to  pipe  diameter: 

P  =  -  (33) 

D 

The  density  is  calculated  using  the  ideal  gas  law  (temperature  is  assumed  constant): 

P 

(34) 


A 

Mo 


u  ' 


The  gas  expansibility  is  determined  as: 


e  =  1 -(0.4 1  +  0.35  */?4)*- 


AP 


K*  P, 


The  discharge  coefficient  for  flange  taps  is  calculated  as: 
C  =  I 


0.598  +  0.468  *  (,tf4  + 10  */?12)  *  j34)  +  (0.87 +  0M*  fi4) 


1M4) 

Ren 


(35) 


(36) 


The  last  term  in  the  discharge  coefficient  calculation  is  small  compared  to  the  other  term 
and  is  neglected  in  the  code. 


%Matlab  code  used  to  create  oxygen  flow  rate  chart 


%Clear  memory 
clear  all; 
clc; 


%Orifice  plate  hole  diameter  (m) 
d=0. 00568706; 

%Pipe  diameter  (m) 

D=0. 0157988; 

%Beta  ratio 
b=d/D; 

%Orifice  and  pipe  diameters  (mA2) 
At=pi/4*dA2; 

Ap=pi/4*DA2; 

%Coefficients 

k=1.29; 

C=(0.598  +0.468*(bA4+10*bA12))*sqrt(l-bA4); 
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%Average  upstream  pressure  (psig) 

P  210: 

%Convert  to  Pa 
pl=P*6894. 757293 17; 

%Temperature  (K) — assumed  value 
T=297; 

%Gas  constant  for  oxygen  (J/kg-K) 
r=260; 

%Standard  pressure  (Pa)  and  temperature  (K) 
pst=101325; 

Tst=289; 

%Density  (kg/mA3) 
rho=pl/r/T; 

i=l; 

dp(i)=0; 
for  1=1:1:101 
%Expansion  coefficient 
e(i)=l-(0.41+0.35*bA4)*dp(i)/k/pl; 

%Flow  rate  (kg/sec) 

qm(i)=e(i)*C*At*sqrt(2*rho*dp(i))/sqrt(FbA4); 

%Other  flow  rates 

qa(i)=qm(i)/rho; 

qs(i)=qa(i)*pl*Tst/pst/T; 

%Pressure  differential  (Pa) 
dp(i+l)=dp(i)+6894. 757293 17; 
end 

%Print  chart 

lprintf('\n\t\t\t\t\t\tDifferential  Pressure  vs.  Flow  Rate\n') 
fprintf('\t  0\t  l\t  2\t  3\t  4\t  5\t  6\t  7\t  8\t  9\t\n') 
for  i=l :  1 0: 1 00 

(print  f('%2.0  l\t%5 ,4fl%5 .4  ft%5 .4  ft%5 .4ft%5 .4ft%5 ,4ft%5 ,4ft%5 .4ft%5 .4ft%5 .4  l\n 
’,i-l,qm(i),... 

qm(  i+ 1 )  ,qm(  i+2) ,  qm(  i+3  ) ,  qm(  i+4) ,  qm(  i+5  ) ,  qm(  i+6) ,  qm(  i+7) ,  qm(  i+8) ,  qm(  i+9)) 
end 

fprintf('  1 00\t%5 .4f\n',qm(  101)) 
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APPENDIX  D:  UNCERTAINTY  CALCULATIONS 


Regression  Rate 

The  regression  rate  equation  was  reduced  to  Eq.  37: 


(37) 


where  r  is  the  regression  rate,  r,  is  the  initial  combustion  port  radius,  %  is  the  burn  time, 
nij  is  the  initial  fuel  grain  mass,  and  in/  is  the  final  fuel  grain  mass.  The  uncertainty  in  the 
regression  rate  is  given  by  Eq.  38: 


CO,  = 


dr 

v3r(. 


0)„ 


+ 


r  dr  * 
—  0). 

Kdh  ‘y 


+ 


dr 

\drrii 


Q)„, 


+ 


f  dr  V 

dm  C°mf 

v  f  J 


(38) 


where  COr  is  the  uncertainty  in  the  regression  rate,  cor  is  the  uncertainty  in  the  initial 


combustion  port  radius,  cot  is  the  uncertainty  in  the  burn  time,  com  is  the  uncertainty  in 

the  initial  fuel  grain  mass,  and  com  is  the  uncertainty  in  the  final  fuel  grain  mass.  The 

uncertainties  are  given  as:  cor  =  ±0.001  m,  cot  =±0.05  s,  0)m  =±0.001  kg,  and 

( y  =  ±0.05/72  , .  Due  to  spallation  of  unburned  mass  from  the  motor,  the  uncertainty  of 

the  final  unburned  fuel  mass  was  conservatively  estimated  as  5%  of  the  final  mass.  The 
other  uncertainties  come  from  the  limitations  of  the  instrumentation. 

Specific  Impulse 

For  ease  in  performing  the  uncertainty  analysis,  the  specific  impulse  can  be  simplified  to 
Eq.  39: 
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(39) 


'7'  go(mf+rhoJb) 

where  Isp  is  the  specific  impulse,  go  is  acceleration  due  to  gravity  at  sea  level  (9.81  m/s2), 
F  is  the  thrust,  and  mox  is  the  oxidizer  mass  flow  rate.  The  uncertainty  in  the  specific 


impulse  is  calculated  by  Eq.  40: 


® isP  = 


/  dlsp 
v  dF 


+ 


rdlsp  V 
— —  co  , 

Kdmf 


+ 


dlsp 

v  dmox 


CO- 


\2  C 

+ 


dlsp 


(40) 


where  coF  is  the  uncertainty  in  the  thrust  and  ojih  is  the  uncertainty  in  the  oxidizer  mass 

flowrate.  The  uncertainties  are  given  as:  COF  =  ±0.0  li7  and  C0lh  =±0.25 mox.  The 

thrust  uncertainty  comes  from  the  instrumentation  and  the  oxidizer  mass  flow  rate 
uncertainty  is  the  average  from  the  uncertainty  due  to  the  pressure  transducers. 
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APPENDIX  E:  INSTRUMENTATION  RECOMMENDATIONS 


Differential  Pressure  Transducer  for  Oxidizer  Mass  Flow  Rate  Measurement 
An  Omega  PX750-100-HDI-SS  differential  pressure  transducer  is  recommended.  At  the 
time  of  this  writing,  this  transducer  costs  $1755.  This  transducer  will  also  require  a 
oxygen  cleaning  service  to  allow  it  to  operate  with  oxygen. 

Thermocouple  for  Exit  Temperature  Measurement 

To  successfully  measure  the  nozzle  exit  plane  temperature,  a  thermocouple  is  needed  that 
can  survive  in  a  high  speed,  high  temperature  (at  least  1800  K),  and  oxidizing 
environment.  A  R  or  S  type  thermocouple  and  a  platinum-rhodium  sheath  are 
recommended.  For  an  exact  price,  a  custom  quote  will  have  to  be  requested  from  Omega. 
Also,  the  response  time  of  the  thermocouple  will  have  to  be  addressed  to  make  sure  the 
burn  time  of  the  rocket  motor  is  sufficiently  long  enough  for  the  thermocouple  to  reach 
saturation. 
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APPENDIX  F:  SUMMARY  OF  TEST  RESUFTS 
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